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CHAPTER  I:  INTRODUCTION 


Wind  generated  or  short  period  waves  continually  arrive  at  the  coast. 

The  approaching  waves  become  unstable  at  a  certain  depth  and  the  tops  of  their 
crests  spill  down  or  plunge  over  their  forward  faces.  The  wave  height 
decreases  as  the  wave  energy  is  converted  into  turbulent  eddies  in  the  surf 
tone.  If  waves  break  at  an  angle  to  the  shore,  they  induce  a  longshore 
current  in  the  surf  zone.  The  current  acts  somewhat  analogous  to  a  river, 
transporting  sediment  mobilized  by  the  breaking  waves.  Coastal  engineers  have 
long  worked  to  correlate  sediment  movement  and  current  velocities  to  predict 
sediment  transport,  shoreline  evolution,  and  pollutant  transport.  This 
''equires  accurate  estimation  procedures,  or  models,  of  the  longshore 
current.  This  report  presents  an  analytical  longshore  current  model  for 
engineering  use.  The  model  employs  an  expression  developed  in  this  report  to 
describe  the  nonlinearity  of  the  wave  height  decay,  and  it  also  includes  the 
effect  of  wave  setup,  finite  incident  wave  angles,  and  lateral  mixing.  The 
advantages  of  an  analytical  model  over  a  numerical  model  are  the  ease  of 
discerning  the  functional  dependencies  of  the  physical  parameters  and  the  ease 
of  applying  the  model. 

Waves  transfer  momentum  from  offshore  to  the  nearshore.  In  the 
nearsnore,  the  waves  break  when  they  reach  a  depth  comparable  to  their  height, 
and  the  wave  energy  is  dissipated  in  the  surf  zone.  Waves  breaking  at  an 
angle  to  the  shoreline  Induce  a  current  parallel  to  the  shoreline  due  to 
changes  in  the  longshore  component  of  momentum.  The  balance  of  momentum  is 
conserved  in  the  surf  zone  by  the  external  forces  of  bottom  and  surface  shea'- 
stresses.  The  change  in  the  onshore  component  of  momentum  also  causes  a 
change  in  the  mean  water  level  in  the  surf  zone,  known  as  wave  setup. 

Momentum  is  also  diffused  or  transported  by  turbulent  eddies. 

Water  motion  in  the  surf  zone  is  extremely  complex.  The  flow  is 
unsteady  and  three-dimensional,  with  dynamic  upper  and  lower  boundaries.  Nc 
adequate  theoretical  description  of  water  motion  in  the  surf  zone  presently 
exists.  Therefore,  to  predict  longshore  currents  it  is  necessary  to  simpliry 
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the  problem  by  considering  an  idealized  environment  and  to  include  a  certain 
amount  of  empiricism.  Applying  various  degrees  of  simplification,  many 
investigators  have  calculated  longshore  currents  analytically  and  numerically 
using  empirical  correlations,  continuity  of  water  mass,  energy  flux,  and 
momentum  flux. 

In  1967  Galvin  reviewed  the  state  of  the  art  of  longshore  current  pre¬ 
diction.  He  concluded  that  the  best  approach  at  the  time  was  the  prediction 
of  longshore  current  velocity  through  empirical  correlation  of  data,  but  he 
cautioned  that  the  available  data  were  not  reliable.  Much  progress  has  been 
made  in  the  prediction  of  longshore  currents  since  the  review  by  Galvin.  The 
progress  mainly  is  due  to  the  introduction  of  the  concept  of  radiation  stress 
by  Gonguet-Higgins  and  Stewart  (1962,  1963,  1964).  Radiation  stress  is  used 
to  calculate  the  flux  of  momentum  parallel  to  the  shoreline  due  to  incident 
waves.  Bowen  (1969),  Longuet-Higgins  (1970a,  1970b),  and  Thornton  (1971)  were 
the  first  to  apply  radiation  stress  concepts  in  the  equations  of  motion  to 
predict  longshore  currents. 

In  the  radiation  stress  approach  it  is  necessary  to  specify  the  wave 
height  through  the  surf  zone  a  priori,  but  the  mechanisms  that  determine  the 
wave  height  in  the  surf  zone  (wave  breaking,  wave  deformation,  and  energy 
dissipation)  are  not  well  'understood.  No  quantitative,  first-principle 
theoretical  model  of  wave  height  decay  exists;  therefore  an  empirical  apprcacn 
13  taken  in  longshore  current  modeling.  The  standard  assumption  is  made  that, 
in  the  surf  zone  (after  initial  wave  breaking),  the  wave  heignt,  H  ,  is 
described  as  a  linear  function  of  water  depth,  h  ,  in  the  form 


H  •  Yh 


where  Y  is  a  constant  of  proportionality.  This  is  known  as  the  spilling 
breaker  assumption  because  it  holds  fairly  well  for  waves  classified  as 
spilling  breakers.  However,  several  investigators  show  that  this  is  not  vaii: 
in  general  (Horikawa  and  Kuo  1967;  Nakamura,  Shirashi,  and  Sasaki  '9c'; 
and  lamfield  1967;  Oivoky,  Le  Mehaute,  and  Lin  1970;  Dally,  Dean,  and 
Dalrymple,  1985a,  1985b),  and  it  is  especially  inappropriate  for  mild  octtim 
slopes,  on  which  waves  tend  to  break  by  plunging. 
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Investigators  following  Bowen.  Longuet-Higgins ,  and  Thornton  built  on 
the  radiation  stress  approach  by  eliminating  some  of  the  simplifying 
assumptions  and  making  more  general  models.  But,  all  have  retained  the 
spilling  breaker  assumption  despite  its  proven  invalidity. 

This  investigation  examines  the  effects  of  using  nonlinear  wave  heignt 
decay,  namely  a  power  law  decay,  on  the  prediction  of  the  longshore  currents 
distribution.  The  power  law  wave  height  decay  is  of  the  form 

H  =  (h/hj^)'’  ;-2: 

where  the  subscript  b  indicates  breaking  conditions,  and  the  exponent  n  ,  t 
be  determined  empirically,  is  assumed  to  be  dependent  on  the  beach  slope  and 
the  breaking  wave  conditions.  It  will  be  shown  in  this  report  that  a  closed- 
form  solution  for  the  longshore  current  distribution  can  still  be  derived  if 
Equation  1-2  is  employed  instead  of  Equation  1-1. 

The  main  body  of  this  report  begins  with  a  review  of  previous  longshore 
current  models.  Special  attention  is  paid  to  the  Longuet-Higgins  model 
because  it  has  served  as  the  basis  for  most  models  that  followed.  Next,  the 
wave  height  decay  portion  of  this  study  is  presented.  Seven  independent  data 
sets  are  empirically  fit  to  the  wave  height  decay  power  law,  and  the  exponent 
of  the  power  law  is  parameterized.  Then,  an  analytical  longshore  current 
model  is  derived  from  the  equations  of  motion  based  on  the  radiation  stress 
approach.  The  effects  of  large  angles  of  wave  incidence  and  of  lateral  mixir. 
a.'"e  included  in  the  model.  The  current  model  gives  the  longshore  current  as 
function  of  distance  offshore,  incident  wave  conditions,  beach  slope,  frictic 
coefficient,  and  a  parameter,  P  ,  expressing  the  relative  importance  of 
lateral  mixing  and  bottom  friction  as  introduced  by  Longuet-Higgins  ;'5'’lb'. 

Review  of  Previous  Models 

The  radiation  stress  approach  to  modeling  longshore  currents  was 
developed  independently  by  Bowen  (1969),  Longuet-Higgins  (19?Ca,  15'^Gd',  and 
Thornton  (1971).  Although  the  three  moaels  are  similar,  there  are  difference 
in  the  ass'umptions  made  in  the  bottom  shear  stress  and  lateral  mixing  te-'me. 


The  former  two  authors  developed  analytical  solutions  for  a  plane  beach;  the 
latter  developed  a  numerical  solution  for  arbitrary  profiles  of  straight, 
parallel  contours,  using  a  more  realistic  bottom  friction  stress.  The 
Longuet-Higgins  moael  is  the  easiest  and  most  straightforward  to  use  (the 
solution  of  Bowen  is  In  terms  of  Bessel  functions  and  the  model  of  Thornton 
requires  a  numerical  solution) ,  and  appears  to  give  very  acceptable  results 
for  a  plane  beach.  The  Longuet-Higgins  model,  therefore,  has  been  used  as  the 
basis  for  more  recent  longsnore  current  models.  A  review  of  the  Longuet- 
Higgins  model  is  given,  followed  by  overviews  of  other  momentum-based 
models.  Basco  (’982)  presents  an  thorough  review  of  surf  zone  current 
literature  with  an  annotated  bibliography  (Basco  and  Coleman  1932).  Table 
gives  an  intercomparison  of  selected  models  of  the  longshore  current 
distribution  across  the  surf  zone. 

Longuet-Higgins .  Longuet-Higgins  (l970a,  1970b)  derives  an  analytical 
model  for  the  steady  longshore  current  from  the  governing  equations  of  water 
motion.  He  makes  the  assumptions  'given  in  Table  1-2;  in  addition  he  assumes 
linear  wave  height  decay  given  by  Equation  1-1.  The  equation  of  motion  for 
the  longshore  direction  for  this  Idealized  case  reduces  to  a  balance  between 
the  local  wave  stress,  the  stress  due  to  horizontal  turbulent  eddies,  and  the 
time-averaged  bottom  friction  stress.  The  local  wave  stress  is  the  driving 
force  of  the  currents,  and  it  is  the  net  stress  in  the  longshore  direction 
exerted  by  the  waves  on  the  water  in  the  surf  zone.  This  stress  is  calculated 
from  the  radiation  stress.  The  bottom  shear  stress  is  linearized  by  assuming 
tne  incident  wave  angle  is  small  and  the  steady  current  is  weak  compared  witn 
the  -wave  orbital  velocities.  These  ass'jinptions  reduce  the  bottom  shear  stress 
to  the  product  of  the  orbital  velocity  and  the  longshore  current  speed.  The 
lateral  mixing  stress  is  a  function  of  the  horizontal  eddy  coefficient. 
Longuet-Higgins  assumes  the  horizontal  eddy  coefficient  is  proportional  to  tne 
offshore  distance  multiplied  by  a  typical  velocity,  the  shallow-water  wa /e 
celerity.  The  distances  (measured  from  the  mean  shoreline)  are 
nondimensional ized  by  the  distance  from  the  mean  shoreline  to  the  breaker 
line.  The  longshore  current  velocity  is  nondimensional ized  by  the  velocity 


•  •  • 

y  %"  V 


k-’w-Vj 


AW 

m 


•.'N 


^  / 

•  n>\' '. 


."J 


10 


V 


TUBLE  l-l 


Table  i-2 


Longshore  Current  Niodel  Ass'jmptiona 


WAVE  FIELD 


Monochromatic  waves 

Linear,  shallow-water  wave  theory 

Steady  state  wave  field 


3  EACH 


Plane,  sloping  beach 

Impermeable  beach 

Hydrostatic  pressure  distribution 


FLUID 


Incompressible,  homogeneous  fluid 


CUaSENT 


Current  constant  through  depth  and  time 
Current  homogeneous  in  the  longshore  direction 
Current  weak  relative  to  the  wave  orbital  velocity 


•VECLECTED  STRESSES 


No  wind  stress 

No  atmospheric  pressure  gradient 
No  wave-current  interaction 
No  Coriolis  force 
No  tide 
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at  the  breaker  line  when  the  effect  of  lateral  mixing  is  omitted.  The  stress 
balance  is  described  by  a  second-order  differential  equation  with  a  closed- 
form  solution.  The  solution  is  a  function  of  the  relative  effects  of  lateral 
mixing  and  bottom  friction.  Longuet-Higgins  does  not  include  wave  setup 
explicitly,  but  he  suggests  modifying  the  beach  slope  to  include  the  change  in 
water  depth  due  to  wave  setup.  He  also  does  not  include  refraction  because 
the  angle  of  wave  incidence  is  assumed  small. 

The  strong  points  of  the  Longuet-Higgins  model  are;  (a)  the  model 
solution  is  simple  and  easy  to  apply  and,  (b)  the  model  results  compare  well 
to  availaole  data.  The  weak  points  of  the  model  are;  (a)  the  numerous 
simplifying  assumptions,  and  (b)  the  spilling  breaker  assumption  in  the 
lateral  mixing  and  bottom  stress  terms  was  applied  seaward  of  the  breaker  line 
where  it  is  no  longer  valid. 

Bowen.  The  Bowen  (1969)  model  differs  from  the  Longuet-Higgins  model  in 
several  ways.  Bowen  assumes  the  bottom  shear  stress  is  proportional  to  the 
longshore  current  speed,  neglecting  the  contribution  of  the  wave  orbital 
velocity.  He  also  does  not  account  for  the  effect  of  variation  in  depth  in 
the  lateral  nixing  stress.  Although  Bowen  simplifies  the  stress  terms 
considerably  more  than  Longuet-Higgins,  his  solution  is  more  complicated.  The 
solution  is  in  terms  of  Bessel  functions  and  is,  therefore,  more  difficult  to 
use.  On  the  positive  side,  Bowen  explicitly  includes  wave  setup  in  the  surf 
cone,  and  he  neglects  it  outside  the  surf  zone  where  it  is  negligible  compa''ed 
to  the  depth. 

Thornton.  Thornton  (1970  uses  solitary  wave  theory  in  the  surf  zone  to 
specify  wave  celerity  and  linear  wave  theory  outside  the  surf  zone.  Thornton 
relaxes  the  plane  beach  assumption,  but  still  assumes  a  beach  of  straight  and 
parallel  contours.  He  also  includes  setup  and  refraction  inside  and  outsi-de 
the  surf  zone.  Thornton  uses  Prandtl's  mixing  length  hypothesis  to  calculate 
the  horizontal  eddy  coefficient  in  the  lateral  mixing  stress.  He  assumes  the 
horizontal  eddy  coefficient  is  equal  to  the  amplitude  of  wave  particle  motion 
multiplied  by  water  particle  velocity  fluctuations  due  to  waves  in  the  snore 
normal  direction.  The  Jonsson  (1967)  friction  factor  for  turbulent  flow  was 
ised  in  the  bottom  stress  term.  Thornton  also  does  not  account  for  the 
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variation  in  depth  in  the  lateral  mixing  stress.  Thornton's  model  requires  a 
numerical  solution. 

James.  James  (197^)  uses  hyperbolic  wave  theory  in  the  surf  zone  and 
linear  wave  theory  far  outside  the  surf  zone  with  a  transition  region  in 
between  to  calculate  the  wave  stress.  Hyperbolic  wave  theory  is  an 
approximation  of  cnoidal  wave  theory  which  is  believed  to  describe  the  wave 
form  in  the  surf  zone  better  than  linear  theory.  James  includes  refraction, 
setup,  and  return  flows  (to  insure  the  mean  shoreward  mass  flux  is  zero).  He 
also  eliminates  the  weak  current  assumption.  Outside  the  surf  zone,  he  uses 
experimental  results  to  define  the  eddy  coefficient  to  be  proportional  to  tne 
inverse  of  the  depth.  James  relaxes  the  plane  beach  assumption,  but  requires 
the  beacn  slope  to  be  mild.  The  mild  slope  assumption  mav  invalidate  the 
linear  wave  height  decay  assumption  (as  stated  earlier).  -,*so,  the  model  is 
formulated  as  a  set  of  differential  equations  that  must  be  solved 
numerically.  This  model  is  much  too  complicated  for  practical  engineering 
use . 

Jonsson,  Skovgaard,  and  Jacobsen.  Jonsson,  Skovgaard,  and  Jacobsen 
return  to  using  linear  wave  theory  throughout  the  nearshore  region. 
They  use  a  nonlinear  bottom  shear  stress  and  introduce  a  friction  factor  tnat 
is  an  interpolation  between  the  friction  factor  for  waves  only  and  tne 
f-'iction  factor  for  currents  only.  Jonsson,  Skovgaard,  and  Jacobsen  adopt 
Thornton's  ()97l)  formulation  for  the  lateral  mixing  stress,  but  they  oo 
account  for  the  variation  in  depth.  The  model  is  a  differential  equation 
whicn  is  solved  numerically. 

Keeley  and  Bowen.  Keeley  and  Bowen  (1977)  take  into  account  longshore 
variations  in  longshore  currents,  removing  the  assumption  of  the  current  bein 
homogeneous  in  the  longshore  direction.  Spatial  variations  in  the  longshore 
current,  typical  in  the  field,  are  caused  by  irregular  bathymetry  and  spatial 
variations  in  the  wave  field.  Keeley  and  Bowen  follow  the  Longuet-H iggi ns 
der’ivation  of  the  longshore  current  due  to  obliquely  incident  waves,  but  om:t 
the  lateral  mixing  stress.  They  linearly  add  the  currents  due  to  obliquely 
.ncident  waves,  variations  in  the  wave  height  in  the  longshore  direction, 
variations  in  the  wave  angle  in  the  longshore  direction,  and  nonlinear  effect 


'due  to  the  advection  term  in  the  longshore  momentum  balance).  They  also 
include  wave  setup.  The  Keeley  and  Bowen  model  must  be  driven  by  a  refraction 
model  which  provides  the  variation  of  wave  heights  and  angles  in  the  lohgshore 
direction.  The  contributions  of  the  longshore  variation  in  wave  height  and 
the  nonlinear  effects  to  the  longshore  current  are  small.  The  model  requires 
a  numerical  solution. 

liu  and  Dalr’/mple.  Liu  and  Dalrymple  (■9’3)  present  a  weaK  current 
model  and  a  strong  current  model.  Both  models  include  the  effects  of  large 
incident  wave  angle  and  wave  setup,  but  exclude  the  lateral  mixing  stress.  In 
the  weaK  cur''ent  model,  the  longshore  current  velocity  is  assumed  small 
compared  to  the  wave  orbital  velocity.  In  the  strong  current  model,  the 
longshore  current  is  assumed  to  be  of  the  same  order  of  magnitude  or  larger 
than  the  wave  orbital  velocity.  The  absolute  value  of  the  total  velocity 
.'longsr.ore  current  plus  wave  orbital  motion)  is  approximated  with  a  truncated 
binomial  series.  For  the  weaK  current  model,  the  bottom  stress  term  is 
simplified  to  a  linear  function  of  the  current  velocity  using  the  weaK  current 
assumption.  The  solution  of  the  weak  current  model  is  in  closed  form.  The 
strong  current  model  results  in  a  nonlinear  ordinary  differential  equation 
solved  numerically.  The  solution  of  the  strong  current  model  is  found 
iteratively  because  the  setup  is  not  known  a  priori.  The  neglect  of  lateral 
mixing  -imits  the  use  of  this  model. 

Kraus  and  Sasaki.  Kraus  and  Sasaki  i'979a,  '!9'’9b)  add  still  anotner 
improvement  to  the  lineage  of  momentum-based  longshore  current  models.  Thei.- 
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model  includes  the  effects  of  large  Incident  wave  angles  and  one  lateral 
mixing  stress  (omitted  by  Liu  and  Dalrymple).  They  assume  that  the  magr.iti.de 
of  the  longshore  current  is  small  compared  to  the  wave  orbital  velocity. 

Setup  is  approximated  by  modifying  the  beach  slope  as  suggested  by  i.orguet- 
Higgins.  Similar  to  the  Liu  and  Dalrymple  strong  current  model,  Kraus  and 
Sasaxi  approximate  the  absolute  value  of  the  total  velocity  'wave  orb.tal  pl..s 
longshore  current)  with  a  truncated  binomial  expansion.  Inside  the  s,rf  cine, 
tney  also  apply  the  approximation 
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derived  from  a  trigonometric  identity,  shallow-water  approximations  for 
wave  celerity,  and  Snell's  law,  where  9  is  the  angle  of  wave  incidence 
model  has  an  analytic  solution  in  the  form  of  an  infinite  series  of 
successively  smaller  terms.  Kraus  and  Sasaki  verified  the  model  with 
laboratory  data  (Mizuguchi  et  al.  1973)  and  thei"  own  field  data. 
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CHAPTER  II:  WAVE  HEIGHT  GECAY 

Th-3  derivation  of  the  wave-induced  longshore  current  requires  knowledge 
of  the  wave  height  and  the  gradient  of  the  wave  heignt  in  the  surf  zone. 
Historically,  the  wave  height  in  the  surf  zone  has  teen  estinated  as  a  linear 
function  of  the  water  depth, 


support  tne  ass'uraption  with  laboratory  data  on  a  slope  of  i/l2.  The  Y-values 
ranged  from  0.9  to  1.3.  This  empirical  expression  is  attractive  because  of 
;t3  simplicity,  but  the  surf  zone  wave  height  decay  is  not  linear  in  general 
as  nas  been  noted,  for  example,  by  Horikawa  and  Kuo  ;'967),  Street  and 
Gamfield  ;'96'''',  and  Van  Dorn  ['■97'''  on  the  basis  of  their  carefully  perfcrme- 
laboratory  experiments.  Figure  2-i  shows  idealized  curves  fit  to  laboratory 
uave  height  decay  data.  The  curves  are  increasingly  concave  upward  with 
:ecrea3:r.g  beacn  slope.  The  purpose  of  this  chapter  os  to  develop  an  ernc.-ioal 
power’  law  decay  model  to  describe  the  wave  heignt  decay  mere  accurately  than 
tne  linear  model,  but  still  retain  the  useful  simple  form  of  tne  linear 
model.  The  simple  form  will  allow  the  longshore  current  mcdel  to  oe  solve: 
analytically. 

The  dissipation  of  wave  energy  in  the  surf  zone  is  d..e  primarily  t: 
turbulence  (Horikawa  and  Kuo  1  967;  Sawaragi  and  Iwata  '?’5;  hio-guoni  'r'O’; 
Gaily,  Dean,  and  Dalrymple  1985a,  i985b;  and  others;.  The  power  law  oeiay 
mcdel  is  entirely  empirical.  The  model  is  not  meant  to  -eplace  more 
sophisticated  models  based  on  the  physics  of  the  turbulent  energy 
dissipation.  These  more  sophisticated  models  solve  tne  energy  fl_x  eouat. in 
in  tne  surf  zone, 
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lere  is  the  energy  flux  and  e  is  tne  e'',e''gy  d.ssioat. 
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For  completeness,  some  of  these  models  are  briefly  described.  Le  riehaute 
;i963)  approximates  a  breaking  wave  as  a  hydraulic  jump,  substituting  the 
energy  dissipation  of  a  hydraulic  jump  for  e  in  Equation  2-2.  The  same 
approach  with  some  variatiohs  is  applied  to  periodic  laboratory  waves  by 
uivoky  et  al .  (1970),  Hwang  and  Divoky  (1971),  and  Svendsen  ;i98'i,  1985). 
Battjes  and  Janssen  (1979)  also  use  the  hydraulic  jump  model,  but  apply  it  to 
random  laboratory  waves.  Thornton  and  Guza  (1983)  refine  the  approach  of 


Battjes  and  Janssen  and  apply  it  to  both  laboratory  and  field  data. 


Although  the  hydraulic  jump  model  appears  to  give  the  best  explanation 
of  the  physios  of  wave  breaking,  three  other  approaches  are  mentioned  because 


of  their  uniqueness  and  insight.  Horikawa  and  Kuo  (1967)  model  surf  zone 
energy  dissipation  due  to  bottom  friction  and  turbulence  using  solitary  wave 
theory.  The  turbulence  is  assumed  to  decay  exponentially  with  distance  from 
the  break  point.  The  results  are  good  for  a  horizontal  bed,  but  poor  for  a 
plane  sloping  bed.  Hizuguchi  (1981)  models  the  surf  zone  energy  dissipation 
by  replacing  the  molecular  viscosity  with  the  turbulent  eddy  viscosity  in  the 
solution  for  internal  energy  dissipation  due  to  viscosity.  Hizuguchi's  model 
allows  more  complex  beach  profiles  (step-type  beaches)  and  reformation  and 
second  breaking  of  waves.  The  model  gives  good  results  when  tested  with 
laboratory  data  for  wave  breaking  on  a  horizontal  beach,  a  1/10  slope  plane 
beach,  and  a  step-type  beach.  But,  Hiruguchi  admits  that  the  eddy  viscosity 
assumption  is  "obscure."  The  model  requires  a  numerical  solution.  Dally  et 
al.  :'935a,  ’985b)  propose  what  they  call  an  intuitive  approach.  The 


iissipation,  c  ,  in  Equation  2-2  is  assumed  to  be  proportional  to  the 


difference  between  the  local  energy  flux,  ECg,  and  the  "stable"  energy  flux, 


(k/h)  (EC,  -  e; 

i 


where  <  is  a  dimensionless  decay  coefficient  and  h  is  tne  local  still¬ 


water  depth.  The  stable  energy  flux  is  found  to  be  associated  with  a  wave 
reight  equal  to  approximately  0.35  to  O.UO  times  the  local  depth.  This 
apprcacr.  allows  a  breaking  wave  to  stabilize  or  reform  and  break  again.  T-e 
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formulation  also  allows  for  an  arbitrary  beach  profile  and  the  inclusion  of 
wave  setup,  but  this  requires  a  numerical  solution.  Analytical  solutions  are 
derived  for  simple  profiles  (horizontal  bottom,  sloping  bottom,  and  Dean's 
;i977)  equilibrium  profile).  Results  are  good  in  comparison  to  laboratory 
data.  Since  this  approach  is  so  successful,  the  power  law  decay  model  will  be 
compared  to  it. 


Power  Law  ?odel  of  Wave  Height  Decay 
In  this  study,  the  wave  height  decay  is  expressed  as  the  power  law 

H  -  Yh^ 

b  b 

This  form  was  chosen  because  it  is  similar  to  the  linear  wave  height  decay 
model,  and  it  reduces  to  the  linear  decay  model  (Equation  2-1)  for  an 
exponent,  n  ,  equal  to  l.O.  Equation  2-4  is  applicable  from  the  breaker  lin 
to  the  mean  shoreline.  Two  constants,  Y  and  n  ,  must  be  specified  in 
Equation  2-^.  It  is  noted  that  the  formulation  of  Dally  et  al.  (1995a,  '935b 
also  requires  specification  of  two  parameters  through  empirical  consideration, 
Tne  importance  of  beach  slope  in  the  decay  profile  is  clearly  shown  in  Figure 
2-'.  Hcrikawa  and  Kuo  also  suggest  the  importance  of  the  wave  steepness, 

1,,  wnere  is  the  deepwater  wave  height  and  is  tne  deepwater 

-.avelength,  and  the  breaking  wave  conditions  (Hj^/h^j)  on  the  decay  profile. 
Fci.lowir.g  a  description  of  the  wave  height  decay  data,  the  procedures  used  tc 
analyze  the  data  and  quantify  Y  and  n  are  explained. 

Seven  independent  sets  of  laboratory  and  prototype  scale  data  ccmprisin 
■35  experimental  runs  on  slopes  of  1 /90  to  WIO  are  used  to  quantify  the  wave 
neight  decay.  These  data  sets  were  obtained  through  a  comprehensive  search  c 
the  literature  in  English  and  Japanese.  Table  2-1  summarizes  the  data.  The 
c''eak..ng  wave  heights  (of  monochromatic  waves)  range  from  4.5?  cm  to  '.3'’ 
and  the  wave  periods  range  from  1.2  3  to  9.0  s.  The  wave  steepnesses 
a-e  between  O.CC31  and  0.091.  The  data  are  listed  in  Appendix  A. 
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Table  2-1 


Data  Summary  for  Wave  Height  Decay 


Source _  Slope  Number  of  R'uns 


Horikawa  and  Kuo  (1967) 
and  Kuo  ('965) 

1  /80 

1  /65 

1  /30 

1  /20 

57 

16 

19 

21 

Maruyama  et  al.  (1983) 

1 /62.5 

, 

1/^5. 5 

1 

1/29.11 

1 

1 /22.2 

1 

Mizuguchi  (1981) 

1/10 

1 

•*  *1^  ' 
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Horikawa  and  Kuo  performed  their  experiment  in  two  parts.  The  1/20  and 
1 /30-3lope  data  were  collected  in  a  flume  17  m  long,  0.7  m  wide,  and  0.6  m 
deep.  The  slope  was  covered  with  a  smooth  rubber  mat.  The  1/55  and  i/So- 
slope  data  were  collected  in  a  flume  75  m  long,  1.0  m  wide,  and  1.2  m  deep. 

The  slope  was  concrete.  The  Maruyama  et  al.  data  were  collected  in  a 
prototype-scale  flume  250  m  long,  S.**  a  wide,  and  1.2  m  deep.  The  initial 
slopes  (1/22.2,  1/29. 'i,  1/^5. 5,  and  ’/62.5(  were  formed  of  sand  and  were, 
therefore,  not  constant  throughout  each  run.  The  Mizuguchi  data  were 
collected  in  a  wave  basin  15  m  long  and  15  m  wide,  but  the  width  was  truncate: 
to  9  m.  The  1 /50-3lope  Saeki  and  Sasaki  data  were  collected  in  a  flume  2^  n 
long,  0.8  m  wide,  and  0.3  m  deep.  The  l/90-3lope  Sasaki  and  Saeki  data  were 
collected  in  a  flume  24  m  long,  0.6  m  wide,  and  l  m  deep.  The  slope  in  both 
cases  was  formed  of  smooth  plastic.  The  Stive  data  were  collected  at  two 
scales  to  compare  scale  effects.  The  large  flume  was  233  long,  5  m  wide, 
and  7  m  deep.  The  slope  was  sand  with  an  initial  slope  of  l/UO.  The  small 
flume  was  55  m  long,  1  m  wide,  and  1  m  deep,  with  a  concrete  slope  of  1/40. 

The  Van  Dorn  data  were  collected  in  a  flume  24  m  long,  0.5  m  wide,  and  at  a 
still-water  depth  of  36  cm.  The  slopes  were  formed  of  plate  glass. 

The  parameter  Y  is  defined  as  the  ratio  of  the  wave  height  to  the 
local  water  depth  at  breaking. 
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by  solving  Equation  2-4  for  Y  with  and  h-h^^.  This  ratio  is  very 

3;gnificant  because  it  specifies  where  a  wave  will  break.  This  is  important 
'.n  the  design  of  coastal  structures,  so  the  specification  of  y  nas 
stimulated  much  Interest. 

McGowan  (1891)  calculates  the  critical  H/h  ratio  for  wave  breaicng  from 
solitary  wave  theory.  His  value. 


This  ratio  is  very 
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gives  a  -easonable  average  of  the  measured  Y-values  from  the  data  S'u-mr.ar  i  oe : 


V- . 


-TV  ^ 


in  Table  2-1,  but  it  does  not  explain  the  measured  variation  in  v  .  Figure 
2-2  shows  McGowan's  expression  versus  measured  values  of  'i'  available  from  tr.e 
data  set  (Appendix  A).  Galvin  (1969)  includes  the  effect  of  beach  slooe,  m  , 
in  his  empirical  relationship 


b  b  b 


where  3,  =  0.92  for  m  >  0.07  and  6.  =  1.90  -  5.85  m  f or  m  <  Q.C”. 

0  b 

and  flier  (1969)  also  include  beach  slope  in  their  empirical  express:: 


Hjj/h^  =  0.72  ♦  5.6  m 


Galvin,  and  Collins  and  Wier  predict  increasing  Y  with  increasing  slope. 

The  data  indicate  that  this  trend  is  correct  (Figures  2-3  and  2-9),  but  the 
large  variations  in  Y  for  a  given  slope  are  not  accounted  for  by  these  two 
equations.  Goda  (1970)  gives  Hjj/h,j  is  terms  of  and  slope  in  graphical 

form  (Figure  2-5)  obtained  by  fitting  field  and  laboratory  data.  Weggel 
(1972)  gives  a  similar  expression  in  terms  of  (where  T  is  tr.e  wave 

period'  and  beach  slope 

H,^/h„  -  b(m)  -  a(m)  H. /T^ 
b  b  b 

-l9m  -19.5m 

where  a.'m)  =  1.36  (1  -  e  )  and  1/b(m)  »  0.59  (1  ♦  e  The  units  :f 

a  ml  are  sec  /ft.  Figure  2-6  shows  Weggel' s  expression  versus  the  measure: 
values  of  Y  .  Singamsettl  and  Wind  (1980)  and  Sunamura  (198')  include  tne 
effects  of  beach  slope  and  wave  steepness  (Hg/Lg)  in  expressions  fc-  n_  h_ . 
oingamsetti  and  Wind's  equation, 

H  /h  •  1.16  (m/(H  /L  )'''2)0-22 
b  b  00 

13  plotted  in  Figure  2-7  against  the  data  summarized  in  Table  2-'.  Sunar.^ns's 
ecuat i on , 


H,/h^  =  1  .1  (m)’''^  (H^/L^)"’ ^ 

b  b  00 
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Figure  2-3.  Comparison  of  Y  calculated  by  the  expression  of 
Galvin  (1969)  and  experimental  results 
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2-u.  Comparison  of  Y  calculated  by  the  expression  of 
(1969)  and  experimental  results 
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rigure  2-1.  '-onparisop.  of  y  calcalated  by  bhe  expression  of 
Singamsebti  and  Wind  (1930)  and  experimental  results 


Figure  2-8.  Comparison  of  y  calculated  by  tr.e  expression 
Sunamura  (1931)  and  experimental  results 


i.3  plotted  In  Figure  2-8  against  the  data.  Equations  2-5  and  2-7  are  quite 
similar.  The  equations  give  a  rough  estimate  of  the  measured  value  of  v 
but  they  do  not  explain  all  of  the  measured  variation  in  v. 

Plots  of  the  measured  Y  versus  (Figure  2-9^  and  beach  slope 

'Figure  2-10)  show  that  Y  decreases  with  increasing  wave  steepness  and 
increases  with  increasing  slope.  Figure  2-il  is  a  plot  of  y  versus  a 
combination  of  these  two  parameters,  known  as  the  surf 

similarity  parameter  (3attjes  i975).  The  plot  shows  some  increase  in  « 
an  increasing  surf  similarity  parameter,  but  the  relationship  is  weak  arc  t 
Cata  are  very  scattered.  Obviously  mucn  effort  has  been  expended  in  the  pa 
to  determine  K,^/h^.  Although  none  of  the  expressior^s  presented  gives  an 
excellent  fit  to  the  data,  the  expressions  of  Singamsetti  and  Wind,  and  of 
Sunaraura  provide  the  best  predictions.  It  is  evident  that  both  the 
measurement  and  the  process  of  wave  breaking  are  very  complex  and  that  the 
phenomenon  has  a  large  variability. 

The  second  parameter  needed  to  quantify  the  wave  height  decay  in  this 
study  is  the  exponent  n  .  The  n-value  in  Equation  2-ii,  obtained  as  a  best 
fit  to  each  of  the  decay  profiles  in  the  data,  was  calculated  by  regression 
analysis.  Equation  2-4  is  nonlinear,  but  it  was  transformed  to  a  linear  fo 
using  natural  logarithms,  and  the  curves  were  fit  to  the  data  by  tr.e  method 
least  squares  (Viller  and  Freund  ’977).  The  method  of  least  squares  minimi 
the  sum  of  the  squares  of  the  vertical  distances  from  the  data  points  to  tn 
-egression  curve.  Equation  2-4  transforms  to 

In  H  -  In  (Yh.  )  *  n  In  (h/h,  ) 
b  b 

in  which  In  is  the  logarithm  to  the  base  e  .  The  previous  equation  .s  ; 
the  form 


Y  -  a  ♦  b  X 


wr.ere  r.  is  unknown.  The  (Yh,  }-term  can  be  treated  either  as  a  <ncwr. 

b 

unknown  in  the  analysis.  The  value  of  vh^  f or  a  particular  run  is 
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equal  to  ,  which  is  available  in  the  data  set.  If  the  value  of  :s 

calculated  by  the  regression,  this  gives  a  better  fit  of  the  curve  to  the 
data.  But,  by  letting  '<  "float"  as  a  free  paraxeter,  it  .s  no  longe-  equal 
to  ,  Therefore,  tne  value  of  ^h^^  is  set  to  for  each  run,  arc 

only  the  value  of  n  is  allowed  to  va""/  to  fit  the  data. 

The  first  attempt  to  quantify  n  was  to  calculate  it  for  each  data  run 
«:tn  tne  regression  analysis,  and  to  plot  these  n-val-.e3  against  significant 
03'’ar.eter3  for  respective  runs.  The  parameters  chosen  to  '•elate  with  the  n- 
val.es  were  beach  slope,  wave  steepness,  surf  si.r.ilarity  parameter,  and  3..rf 
sin.larity  cara.T.eter  at  wave  creaking,  ,  Fig'ure  2-’ 2  shows 

/ersus  beach  slope  for  eucr.  '•un.  The  plot  shows  an  inverse  or  hyperbolic 
-el ationsnip  between  n  and  slope,  but  there  is  much  scatter.  An  n-value 
greater  than  i.O  indicates  a  concave  upward  decay  profile,  so  larger  n-values 
associated  with  smaller  slopes  fits  the  trend  in  Figure  2-1.  Fig’ure  2-13  is  ^ 
plct  of  n  versus  No  correlation  between  n  and  H^/Lq  is  obvious. 

Fig'^^res  2-’*^  and  2-'’ 5  are  plots  of  n  versus  the  deepwater  surf  similarity 
parameter  and  the  surf  similarity  parameter  at  breaking.  The  correlation 
takes  a  hyperbolic  shape  in  both  cases.  These  four  plots  show  that  n  and 

do  not  have  a  strong  relationsnip,  but  n  and  slope  are  rei-ated,  as  was 

V  J 

<ncwn  f-om  the  onset.  The  plots  show  that  for  steep  bottom  slopes,  the 
'.-value  is  lower  and  less  variable,  whereas  on  gentle  bottom  slopes,  tne  n- 
value  is  extremely  variable. 

This  first  attempt  was  encouraging,  but  not  cohcl.^3ive.  The  n-val..e3 
for  the  individual  runs  on  small  slopes  were  extremely  variable.  Iloser 
examination  of  these  runs  (Horikawa  and  Kuo  1/30  and  1  ^65  slopes)  shewed  tnat 
some  runs  had  as  few  as  four  data  points,  and  in  some  cases  the  data  scanned 
only  one-third  of  the  surf  zone  (from  the  breaker  line  i.nshore..  In  'ur.s  witn 
few  data  points,  the  n-values  were  higher.  The  lack  of  inshore  data  points  . c 
the  decay  profiles  evidently  biased  the  re.sults. 

To  eliminate  the  problem  of  sparse  data  in  some  runs,  all  data  f.'  eacn 
slice  were  nondimensionalized  and  combined.  The  wave  height  was 


yyy 

3  to  *  4  '  ■ 


.•  irf* 

yyy 

•■.•Vv 


< 

A'A 

■Jvv 


:gure  2-'i.  Power  law  exponent  versus  wave  steepness 
xperimental  results  'exponent  ieterr.ineo  from  eacr. 
xper i.-nental  run' 


nondimensional Ized  by  the  breaking  wave  height  and  the  water  depth  was 
nondlmenaional ized  by  the  water  depth  at  breaking.  Equation  2-‘i  becomes 

H/Hjj  -  (h/h^)^ 

The  equation  was  transformed  to  the  linear  form 


In  •  n  In  'h/h^' 


and  a  •'egression  analysis  was  performed.  The  results  of  the  regression  a-e 
given  .n  Table  2-2.  The  n-value  is  the  result  of  the  regression  for  tne 
l-.nped  nondimensional  data.  The  values  of  y  and  are  tne  averages  of 

tne  comtined  -uns.  Figure  2-i6  is  a  plot  of  tne  "lumped"  n  versus  the  inverse 
of  bottom  slope.  A  linear  regression  of  n  as  a  function  of  slope  gave  the 
equation 

n  -  1 . 35  ♦  0.009/m 

w.th  a  oor-eiation  coefficient,  r  ,  equal  to  O.^i.  a  regression  was  also 
done  on  a  subset  of  the  data,  the  Horikawa  and  Kuo  data,  and  gave  the  equation 

n  -  0.39  *  O.OiT.'m 
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...  tn  a  oorrelation  coefficient  equal  to  0.86.  The  Horikawa  and  Kuc  data  set 
was  selected  as  a  subset  because  it  Includes  more  runs  and  data  points  per 
slope  than  the  other  data  sets  and  thus  would  be  more  statistically  staole. 
-.gure  2-17  is  a  plot  of  the  "I'umped"  n  versus  the  averaged  L_.  The  plot 
snows  no  obvious  correlation. 

Altnough  the  relationship  between  n  and  beach  slope  is  fairly  clear 
for  the  combined  data,  the  scatter  in  the  data  implies  there  may  be  anctn.er 
.moo-tant  factor,  assuming  that  random  experimental  variability  is  net  tne 
-a.:.''  cause.  To  resolve  the  variation  in  n  for  a  given  slope,  n  versus 

and  n  ve’-sus  ^  were  plotted  for  selected  slopes.  Figures  2-‘;  an.o 
2--}  are  clots  of  n  versus  and  n  versus  y  for  a  slope  of  •21. 

Tne  values  of  n  ,  and  y  are  from  each  of  tne  individual  -un.s.  •'.ot 
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TABLE  2-2 
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averages.  The  n-value  still  shows  no  relation  to  H^/Lq,  tut  n  does  increase 
as  Y  increases.  To  summarize,  n  is  a  function  of  slope  and  y  ,  tut  n 
does  not  appear  to  be  related  to  or  to  the  surf  similarity  parameter  'a 

f unction  of  Hq/L^) . 

As  a  next  step,  to  avoid  the  problem  of  sparse  data  in  some  runs  without 
combining  all  the  data  of  the  same  slope,  runs  with  similar  Y  were  combined 
'so  the  effect  of  Y  is  seen).  Slopes  with  multiple  runs  were  divided  into 
g'-oups  with  Y  centered  on  whole  tenths  of  y  (0.6,  0.7,  0.8,  etc...). 
P.egressions  were  run  on  each  of  the  groups  to  calculate  a  best  fit  n  .  These 
n- values,  the  average  Y  ,  and  the  slope  for  each  group  were  used  as  input  to 
a  multiple  regression  of  n  in  terms  of  Y  and  slope.  The  assumed  form  cf 
tne  equation  is 

n  «  b  ♦  b, /m  +  b.Y  +  b.Y/m 
o  1  2  3 

The  method  of  least  squares  was  used  again  to  fit  a  family  of  curves  to 
the  data.  The  multiple  regression  was  run  on  a  subset  of  the  data,  the 
Horifcawa  and  Kuo  data,  because  it  included  more  data  points  per  slope.  The 
'•esult  of  the  multiple  regression  is 

n  -  0.657Y  *  0.0U3T/m  -  0.0096/m  ♦  0.032  ,2-? 

Equation  2-3  is  plotted  in  Figure  2-20.  The  plot  shows  that  n  :ncrea3es 
witn  increasing  Y  and  decreasing  slope.  The  interaction  term  of  Y  and 
slope  accounts  for  the  increased  steepness  of  the  curves  as  the  slope 
decreases.  Figure  2-21  shows  all  the  data  plotted  against  Equation  2-5.  Tne 
second  term  in  Equation  2-8  is  the  leading  term  for  moderate-to-mild  teaon 
slopes.  The  value  of  n  is  mainly  controlled  by  Y  ,  but  is  also  sensitive 
to  m  because  the  beach  slope  varies  over  an  order  of  magnitude  whereas  the 
value  of  Y  deviates  little  from  unity. 

Figures  2-22  through  2-27  give  examples  of  the  fit  of  the  power  law  wsve 
ne.grt  decay  model  to  the  Horikawa  and  Kuo  laboratory  data.  The  solid  lire 
''ep''esent3  the  power  law  model.  The  n-values  used  were  calculated  f''or. 
Eduation  2-8  with  the  measured  values  of  n  and  beach  slope.  The  dashed  l.re 
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Fig-jre  2-20.  Power  law  exponent  versus  Y  f-om  Equation  2-3  and 
experimental  results  from  Horlkawa  and  Kuo  ^;967)  (exponent 
determined  from  compined  data  with  similar  and  slo-e) 


rigure  2-21.  Power  law  exponent  versus  Y  from  Equation  2-3  and 
experimental  results  (exponent  determined  from  combined  data 
with  similar  Y  and  slope) 
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Figure  2-22.  Comparison  of  wave  height  from  experimental  results,  power 
law  model  (n-1.10,  Y  -l.iS,  and  3lope-l/20),  and  Dally  et  al .  (’985a, 
’935b)  model 
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."igure  2-23.  Comparison  of  wave  height  decay  from  experimental  resu) 
power  law  model  (n-l.6l,  Y  »i.i8,  and  slope-i /20) ,  and  Dally  et  al. 
555a , 1 985b )  model 
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2-2U.  Comparison  of  wave  height  decay  from  experimental  "esulta, 
aw  model  (n-l.uo,  Y  -0.86,  and  3lope-l/30),  and  Oally  et  al . 
’9956)  model 
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gure  2-25.  Comparison  of  wave  height  decay  from  experimental  ■'escl' 
wer  law  model  (n-2.09,  Y  -0.78,  and  3lope-1/65),  and  Dally  et  al . 
995a, 1935b)  model 
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'igure  2-25.  Cor.pariacn  of  wave  height  decay  from  experimental  r 
iz^er  law  model  (n-2, 39.  y  -O-St,  and  slope-1 /65/,  and  Oally  et 
:  1 935a, 1935b)  model 
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Lgure  2-2”.  Comparison  of  wave  height  decay  from  experimetnal  -es 
:wer  law  model  (n-i.75,  ^  »0.6l,  and  slope-i  ^80i  ,  and  Ca*../  et  a* 
'  ?t53 , ' 935b }  model 


or  k  »  0.15  and  the  stable  wave  height  of  Hg  -  3.  *40  h  were  ..sed  td  ca::'.;la 
the  Dally  et  al.  wave  height  decay.  The  decay  profiles  on  the  steep  beach 
slopes  show  almost  linear  decay  with  a  pronounced  setup  '.not  included  in  tne 
Dally  et  al .  curve).  The  decay  profiles  on  the  mild  slopes  snow  a  concave 
upward  shape  as  predicted  by  both  models.  The  Dally  et  al.  model  appears  tc 
fit  the  profile  shape  better  on  the  mild  slopes,  characterizing  the 
'’eformation  of  the  wave.  Overall,  the  power  law  model  gives  a  good  p'-edicti 
of  the  wave  height  decay. 

Tne  power  law  model  predicts  the  wave  height  decay  better  than  the 
linear  decay  model  assumed  in  previous  longshore  current  models.  The  power 
law  model  also  compares  favorably  to  the  more  complex  model  of  Dally  et  al. 
'■995a,  1985b).  The  specification  of  two  parameters,  Y  and  n  ,  is  requir 
:r.  the  power  law  model.  The  parameter  Y  is  the  ratio  of  wave  height  to 
water  depth  at  wave  breaking.  This  parameter  is  best  estimated  by  the 
expression  of  Singamsetti  and  Wind  (Equation  2-6)  or  the  expression  of 
Sunamura  (Equation  2-7).  The  exponent,  n  ,  of  the  power  law  is  3  function 
and  the  beach  slope  as  specified  in  Equation  2-B.  A  closed-form  solutio 
for  the  longshore  current  distribution  is  derived  in  the  next  chapter  using 
the  power  law  model  of  wave  height  decay. 


:HAPTER  III:  DERIVATION  OF  THE  LONGSHORE  CURRENT  DISTRISUTICN  HOCCL 


This  chapter  describes  the  derivation  of  the  closed-form  mathematical 
model  of  the  longshore  current  distribution  based  on  the  power  law  of  wave 
height  decay  in  the  surf  zone  developed  in  Chapter  II,  The  model  is  inter.de 
to  be  ar.  engineering  tod  for  predicting  longshore  cur'rents  and  for  studying 
"el ationsnips  between  physical  factors  generating  the  currents.  The  momentu 
oalance  in  the  longshore  direction  is  the  basis  for  the  model,  but  many 
simplifying  ass'umptions  are  made  in  order  to  provide  a  solution  in  a  form  fo 
practical  use.  The  model  may  be  viewed  as  an  extension  of  the  model  of 
Longuet-Higgins  (1970a,  1970b).  The  effect  of  incident  wave  angles  is 
included  in  the  form  presented  by  Kraus  and  Sasaki  (1978a,  1973b),  but 
truncated  at  second  order  to  allow  easier  application.  The  longshore  curren 
model  is  compared  to  laboratory  data  and  to  the  Longuet-Higgins  model. 

Assumptions 

The  assumptions  used  in  the  derivation  of  the  longshore  curreht  model 
are  listed  in  Table  1-2.  Tnese  assumptions  simplify  tne  mathematical 
development ,  so  an  analytical  solution  becomes  possible.  Similar  assumption 
have  seen  made  in  most  previous  longshore  current  models,  including  numerica 
-coels.  The  assumptions  picture  a  highly  oversimplified  environment,  somew-. 
"emcved  from  the  real  world.  The  longshore  current  is  never  completely  stea 
'see,  e.g..  Meadows  1977),  in  contrast  to  the  steady  state  assumption.  T-.e 
longshore  current  varies  significantly  over  time  periods  as  short  as  c 
minutes.  The  longshore  current  is  also  assumed  to  be  homogeneous  in  t.he 
iongsnore  direction.  Harris  (1969)  describes  this  as  an  "alongshore  system" 
as  opposed  to  a  cellular  system  with  the  longshore  current  feeding  "ip 
currents.  Harris  notes  that  alongshore  s  "’terns  occurred  in  only  'i  percent 
nis  field  observations  performed  on  the  Natal  coast  of  Soutn  Africa,  in 
iac':''ato'’y  wave  basins,  more  conditions  can  be  controlled  e.g.,  the  wave 
field  and  the  beach  slope),  but  the  lateral  boundary  condition  :f  a  homoger.c 


current  in  the  longshore  direction  is  difficult  to  achieve.  If  proper  care  .a 
not  taken,  a  large  circulation  cell  will  tend  to  form  in  a  wave  basin. 

Although  the  assumptions  made  are  restrictive,  the  trends  observed  in  the 
model  are  expected  to  be  applicable  to  more  complex  situations.  If  a  cellular 
system  is  present,  a  circulation  model  (e.g.,  Keeley  and  Bowen  '-977)  should  be 
used.  The  model  presented  in  tnis  report  predicts  the  depth-averaged 
longshore  current  distribution.  This  level  of  sopnistication  is  consistent 
with  most  available  measurements  of  the  longshore  current. 

Equations  of  Motion 

The  equations  of  motion  are  statements  of  Newton's  Second  Law,  conser¬ 
vation  of  momentum.  The  equations  of  motion  for  the  depth-averaged,  steady 
flow  that  must  be  satisfied  are 

u  9u/3x  +  V  3u/9y  -  f^v  =•  -g  Ch*n)  3n/9x  ♦  Il/p)  Z  stresses^  (S"'; 

for  the  x-direction  (shore  normal)  and 

u  3v/3x  ♦  V  3v/3y  +  f  u  =  -g  (h*n)  3ri/3y  ♦  (1/*)  E  stresses  .2-2 

d  y 

for  the  y-directioh  (shore  parallel),  where  u  is  the  mean  current  speed  in. 
tne  x-direotion,  v  is  the  mean  current  speed  in  the  y-direction,  f.,  is  t^.e 
Io''iclis  parameter,  g  is  the  gravitational  acceleration,  and  r  is  tne  mean 
setup  !.-igure  S-M.  A  derivation  of  the  equations  of  motion  is  found  in  Lean 
and  Lalrymple  (198U).  The  continuity  equation, 

"'/3x  [u(n  ♦  h'.Z  *  3/3y  [v(ri  ♦  h)]  -  0 

expressing  conservation  of  mass,  must  also  be  satisfied.  All  y-de’'i  vati  ves 
are  sere  because  of  the  assamption  of  homogeneity  in  the  /-direction  and, 
applying  the  continuity  equation  (Equation  3~3)t  u-0.  The  loriolis  force  is 


v> 


'.egiected.  Equations  S"!  ^nd  3-2  simplify  to 


0  -  -g  (h*n)  3n/9x  ♦  ('/p'  I  stresses 


0  *  Z  stresses 


Equat:cr.  3“*^  '  x-mcment  urn )  is  used  in  Chapter  IV  to  derive  the  wave  setut. 


-quatitn 


is  expanded  to  derive  tne  longshore  current  distribution. 


"ne  stresses  referred  to  in  Equation  3“5  are  the  local  wave  stress, 


,  tne  wind  stress, 


,  the  stress  due  to  lateral  mixing,  3.  ,  and  tne 


f-'idtional  stress  cn  the  bottom,  <By>  .  Equation  3-5  becomes 


y  ■''y 


3,  -  <3> 

'•  y 


The  wind  stress  is  not  included  in  the  general  derivation,  but  it  is  discussed 
later  is  this  chapter. 


local  Wave  Stress 


'he  local  wave  stress  is  the  longshore  force  exerted  on  the  nea.''3nd''e 
water  mass  by  the  incoming  waves,  and  it  is  typically  assumed  to  be  tne  d'l/ 
driving  force  of  the  longshore  current.  The  local  wave  stress  is  cal  dilate: 
^si.ng  the  concept  of  radiation  stress  developed  by  Longuet-H. i ggi  ns  and  Ste-ar 
■'9t2,  ’963,  196^;.  Longuet-Higgi.ns  and  Stewart  define  radiation  st'-ess  as 
tne  excess  flow  of  momentum  due  to  the  presence  of  waves.  The  ri..x  :f  y- 
mcrnentum  parallel  to  the  shore  across  a  plane  x-constant  .3 


F  'sine/O 

X 


'epresents  the  radiation  stress  component  whic.n  is  tne  excess 
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of  x-directed  momentum  in  the  /-direction.  i3  the  energy  flux  in  tne  x- 

iirection  per  unit  distance  alongshore,  C  is  the  wave  celerity,  and  9  is 
the  local  wave  angle.  3y  Snell's  law,  'sine/CI  is  a  constant  and  is 
tnerefore  equal  to  the  same  ratio  at  breaking.  Equation  3-''  can  then  be 
written  as 


S  »  F  isine. /C,  ) 
xy  X  b  b 


Applying  linear  wave  theory,  the  x-directed  energy  flux  is 


F  -EC  cos9 
X  g 


wne^'e  is  the  local  group  celerity  of  the  waves,  and  E 

energy  tensity  per  unit  surface  area.  The  energy  density  is 


is  the  local 


E  -  1/3  p  g  H 

If  waves  do  not  lose  energy  (by  wave  breaking  or  bottom  frictior.:, 
ehe>'gy  flux  is  constant,  but  in  the  surf  zone,  wave  energy  is  certainly 
list.  The  wave  energy  decays  through  the  surf  zone,  and  it  is  zero  at 
3pp-cxiT.ately  the  shoreline.  The  rate  of  energy  dissipation,  D  ,  is 


3F  /9x  -  -D 

X 


The  net  stress  per  unit  area  exerted  by  the  waves  on  tne  water  in  tne 
s^rf  zone  is 


T  »  -  3S  /9x 
y  xy 


J- 

A  '  >  •  •  IT. 


and  from  Equations  3-",  and  3"12, 


r  -  -  s/9x  F  '  sine,.' 
y  X  b 


C  sine^/-^ 


.he  local  wave  stress  is  proportional  to  the  rate  of  energy  dissipation. 
Therefore,  outside  the  surf  zone,  where  energy  loss  is  small  (minimal  wave 
dreaking:  weak  bottom  orbital  velocities,  producing  little  bottom  friction) 
the  wave  stress  is  considered  zero.  I.nside  the  surf  zone,  energy  loss  by  wave 
creaking  is  dominant,  and  bottom  friction  is  also  believed  to  be  significant. 

To  this  point,  the  derivation  is  not  original,  but  has  followed  that  :f 
■_ :nguet-H iggi ns  !'9h0a,  ’9'’Cb).  Next,  the  wave  height  decay  power  law  is 
. nccrpc'"ated  by  applying  Equation  2-^  to  describe  the  broken  wave  heignt  in 
tne  surf  zone.  The  local  wave  energy  ^Equation  3"’u)  becomes 

E  =  '  .'3  p  g  /’'h^.'h'h^.  3-'  - 

and  tne  energy  flux  'Equation  3”9;  becomes 

-  ,  .  n .  ’ 

F  *  1  '3  p  g  Yh.  ih/h.  ,  V  cos 8  '3-‘o 

X  -bb  g 

In  snallow  water,  the  wave  group  celerity  and  the  wave  celerity  are  eqcal  and 
exp-essed  as 

^2  ... 

•g  •  -  ■  * 

by  linear  shallow-water  wave  theory.  Using  Equation  3"'6,  Equation  3"'o 
'’e  d ..  ces  to 


F 

X 


1  /3  p  g 


3/2^2  ;2n»i/2;,,  .'2n-2' 


C038 


wave  setup  is  accounted  for  inside  the  surf  zone  by  altering  tne  oeacn  slice, 
tans  ,  as  suggested  by  Longuet-Higgins  :’9’0a) 

-  dh/dx  •  tanB  •  tan3  '  ( )  *-3^'Sy‘' )  J-'- 


-  me  a' 


'unct  1C 


;f  water  depth  c 


.  ne  «ave  se 


13  assumed  to  be  a 


a 


The  derivative  of  the  energy  flux  in  the  x-direction  is  computed  as  an 
intermediate  step  to  calculating  the  local  wave  stress.  From  Equation  S-'"’ 
and  Equation  3~18 


3F  /9x  -  -tanB  9F  /9h 
X  x 


*  P  'Pn-P'j 

•tans  1  .'3  D  g-  Y  /(h^)  a/3h;h'~“  “  ccse' 


The  Iccal  driving  wave  stress  inside  the  surf  zone  follows  directly  from 
Fixation  3-’9  substituted  in  to  Equation  S'lB,  and  the  wave  stress  outside  t.-.e 


zone,  where  D  is  negligible,  is  zero 


»  ■3/2  2  !  2^-2)  '2''*^/?) 

tana  ’ /9  p  g^  Y‘^/(h^)'  '  9/3h(h'  '  cose)  x  <  x, 

b  b 


(sine^/C^) 


'  3-20 ' 


X  >  x^ 


Lateral  Mixing  Stress 


The  lateral  mixing  stress  is  the  exchange  of  nomentun  caused  by 
nc-izor.tal  turbulent  eddies.  A  review  of  the  influence  of  lateral  mixing  in 
longshore  cur’’ent  modeling  was  recently  made  by  McDougal  and  Hudspeth 


■?3t'.  N’eglect  of  lateral  mixing  predicts  an  unrealistic  discchtinuity  :  r. 


t-e  longshore  current  profile  at  the  breaker  line.  The  lateral  mixing  st-'esa 
emclo/ed  is  of  the  form  used  by  Longuet-Higgins  ;'972b; 


B,  -  3/3x  (pe^^h  3v/3x) 


-here  £.  is  the  lateral  viscosity  coefficient  defined  as  tne  product  of  a 


'•epresentat 1 ve  mixing  length  and  velocity.  The  lateral  viscosity  cceffioient 
..set  is  from  Madsen  et  al .  ii978).  The  representative  mixing  length  -sed  is 


the  distance  to  the  mean  shoreline,  and  the  representative  velocity  is  tne 


-axi'ium  cntital  velocity,  'J 


r  X  U 
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nere  r  13  a  constant.  ihe  maximum  orbital  ve-ocity  is 


U  .  H  C/(2h) 

0 

max 

•  -'2 

*  H  g  "  2  n 

V  l.nea-'  snal  1  ow-water  wave  theory.  Applying  -quation  3“23  arc  tre 
xcressiin  for  the  broken  wave  neight  in  the  surf  cone,  Equation  3-2; 


fx  g  /2  Yh 


,(n-l/2),,,,  .n-1 


he  expression  for  the  lateral  mixing  stress  in  the  surf  zone  from  Equa 
'  anc  3~2u  becomes 


3.  =/5/3x  ;o.'-  X  g'  '~/2  Y  ^  5v/5x:  x  <  x,_ 

ji/jx  iof  X  g  H  h^ 3v/3x)  y  >  x^ 


Bottom  Friction  Stress 


The  bottom  friction  stress  resists  the  flow  along  the  bottom. 


•torn  friction  stress  is  described  by 


B  »  c-  p 


hr 


here,  :<■  is  the  friction  coefficient  and  "J  is  the  total  velocity, 
o-.posec  of  the  wave  orbital  velocity,  1)'^^  ,  and  the  longshore  our-'er: 
eiooity.  The  total  velocity  for  a  longshore  current  system  is 


U  -  LIq  ♦  O.v) 


'3-2'’. 


where  the  notation  (0,  v)  is  used  to  denote  the  x  and  y-oomponents  of  tne 
steady  current.  The  arrows  indicate  vector  quantities.  The  absolute  value  cl 


the  bottom  velocity  is 


ful  =  (IT-  *  2  IT  (0. 


V)  *  v2)’/2 


j--;3  ■ 


Applying  the  assumption  that  v  is  much  smaller  than  'J  and  expanding 
Equation  2-23  with  a  truncated  binomial  series  (retaining  only  first  order 
terms ) ,  yields 


U|  =  |u;|  *  V/  iuj 


( 3-2'; 


fhe  y-component  of  |U|  U  is 


;  IlTi  iT)  -  1"^  I  *  3ln8  v)/lirn  sins  ♦  v) 

'•'o'  n  'n'  a 


(U  I  u  sinS  ♦  V 
'o'  0 


♦  sin^9  v/IlT 
o  '  c 


•y* 

■  *•“  • 


V  U  sine/  U 
o  'o' 


li.nce  the  time  average  of  the  bottom  friction  is  required  to  compute  the  mear 
.ongshore  current,  linear  terms  of  Uq  do  not  contribute  and  can  be  droppeo, 


( iUlU) 


I  V  ( 1  *  Sin's ; 


'3-; 


"ne  resulting  y-component  of  the  bottom  stress  is 


5  -  c  ,0  liT  I V  ( 1  sin  0 )  “^ 


'-■‘I'-V'' . 
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The  time  average  of  the  orbital  velocity  is 


<  U^>  -  (2/n)  U 
0  0 


(3-33 


max 


where  U  is  given  in  Equation  3-23.  Equation  3-32  simplifies  to 


max 


<3  >  -  0  /it  p  V 


/2^(n-1/2; ,  (n-1 )  ,2 


Y  g  h 


H/ 


/(n^) 

b 


'  J  sine) 


1  ♦  3  i  n  9  ) 


X  <  X, 


X  >  X, 


Longshore  Current  Velocity 


Special  case:  small  incident  wave  angle.  Further  simplifications  car. 
oe  made  by  assuming  the  angle  of  wave  incidence  is  small.  The  small-angle 
assumption  also  facilitates  comparison  between  the  present  model  and  the  model 
of  Longuet-Higgins ,  since  Longuet-Higgins  assumes  the  wave  angle  is  small. 

For  this  special  case,  cose  is  approximately  equal  to  unity  and  si.n9 
-.3  approximately  equal  to  zero.  Substituting  the  longshore  stresses 
applicable  to  the  surf  zone  (Equations  3-20,  3-25,  and  3-3^^)  irito  the  stress 
calar.ce  (Equation  3-'3)  gives 


.  *  .  3/2^2  AZx\-\/2)  ,(2n-2)  .  „ 

ta.nS  i4n+1)^l6  p  g  Y  h  /(h.)  sin9^/^, 

b  b  b 


♦  3/3x  (pTYx  g^^^/2  ^3v/3x) 


■j- 


,  ^  1 /2^(n-1 /2) ,n-l 

-  Cj,/Tr  p  V  Yg  h 


Applying  the  plane  beach  assumption  (h  >  tanB  x)  ,  and  simplifying,  Equati 
3-35  becomes 
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'■.  ‘^-  •"•  r^i.'TAr  vr  v- 


r/2  tans*  a/3x  3v/3x)  -  l/ir  c^.  v  x 

NondiJiensionallzlng  x  ,  letting  X  »  x/x^j  ,  results  in 

-  tan8*('4n+i  )/l6  g’'^Y(h  ^(2n-l/2) 

0  D 


(n-1 /2) 


r/2  tans*  3/3X 


' n- ’ /2  ^ 

V/3X)  -  1  /  it  V  X 


(3-3') 


tX  '  XlfjTl 


N" 


If  t'.e  lateral  nfiixing  terii  is  neglected,  the  first  term  to  the  right  of  tne 
eq^al  sign  :n  E:;uation  3-31'  is  zero,  and  the  longshore  current  speed  is  solve, 
for  oi'-ectly 


.-"I--.'-] 


*  1/2  1/2  n 

V  -  tt/c^  tans  (“^n+n/lS  g  (h  )  v  X  sinS^ 


(3-33' 


The  velocity  at  the  breaKer  line  (X  »  l )  for  n  •  1,  neglecting  lateral  mixing, 
IS  defined  as 


5it  tans  s)/2  ,  . 

V  •  -fr  — -  Y(gh.  3in0. 

0  Id  c-  “  b  b 


(3-39' 


following  Longuet-Higgins  (■1970a).  The  value  of  v^  is  the  maximum  possible 
current  speed  for  n  »  i.  As  n  varies,  the  maximum  current  speed  is  given  by 


%ax  -  '^o 


Nondimer.s ional i zing  the  current  speed  v  by  ,  'J  ^  ^ 

Equation  3-37  to 

-(Un-*-1)/5  •  P  a/3X  (X'"*^'^^^3V/3X)  -  V 


■  V  V  w 

n"'*.*  "r 
vw.  **->' 

»•  V  "v  V 


where 


•:-:a 


P  -  (rTr)/(2  c^)  tans 


o  -r  ^ 


Tne  parameter  P  is  nondimens ional ,  and  it  expresses  the  relative  importance 
of  lateral  mixing  (r)  and  bottom  friction  (cf)  • 


54 


.•*v"*\VV 


.  # 


■*'  •  ■ 

>V‘ 


**•'*'  /“  .*•■>•  •  -  •  ■  •  •  '.n  ■*  •  *  -  O  O  S.*  Si*  S  S  S* -s  V*  v*  A.*. 


Calculating  the  derivative  in  Equation  and  rearranging  the  ter'ns 

gives 

P  V”  ♦  ?  {n+3/2)  X  V  -  V  -  -{Un+D/S  X^  rS-j] 

where  the  primes  denote  derivatives  with  respect  to  X,  Equation  3-^3  is  a 
nor.homogeneous  second-order  differential  equation  solved  by  the  method  of 
variation  of  parameters.  The  solution  to  the  differential  equation  is 

V  =  3  xP  *  A  X^  :3--- 

where 

A  =  (4n*l )/5  [1/(1  -  n  P  (2n+1 /2) ) [ 

p  =  -(2n*l)/H  ♦  [( (2n*1  )^/l6)  ♦  1/P]^'’^ 

Substituting  the  longshore  stresses  applicable  outside  the  surf  cone 
'.Equations  3-20,  3-25,  and  3"3^)  into  the  stress  balance  equation  'Equation 
•-o'  yields 

0  -  3/8x  (pr  X  H  8v/3x)  3-^5 

-  Cj./irp  V  g^  (1*3in^9) 

Outside  the  surf  zone  the  wave  height  is  approximated  by  linear  shal 1 ow-wate- 
wave  theory  (Green's  Law)  as 

H  -  (  cose/cose.  )' '^^(h  /h) ' H  3--0 

b  b  b 

Applying  the  small  angle  assumption  and  noting  =■  Yh^  ,  Equation  3'-5 
3..tplifies  to 


✓  y 


Applying  the  small  angle  assumption,  the  plane  teach  assumption,  and  Equation 
3-^'’,  Equation  3~^5  simplifies  to 


5/H  ■ 

0  =  Q  3/5x  x  9v/3x;  -  v/x' 


;  3--8 


r.  ncndimens ional  form 


0  =  Q  3/3X  :X^''‘'3V/3X)  - 


;  3-49 


Q  -  rw/(2  c^)  tans 


■'".e  parameter  Q  is  used  instead  of  ?  seaward  of  the  breaker  line  because 
:r.e  effect  of  wave  setup  is  negligible  in  this  region.  Calculating  the 
lerivative  in  Equation  3”‘^9  and  rearranging  the  terms  gives 


X^V  '  *  5/y  XV'  -  V/Q  =  0 


Eq^3t;:n  3-5O  is  a  homogeneous  second-order  differential  equation  with  tne 
solution 


-1/8  -  (1/514 


Che  quantities  B  from  Equation  3”i4“  and  C  f-om  Equation  3-0'  were 
ootained  by  equating  the  current  and  the  derivative  of  tne  current  inside  ar 
outside  tne  sarf  zone  at  the  breaker  line.  The  general  solution  for  tne 
::ng3no'’e  current  distribution  assuming  a  small  incident  wave  angle  is 


*  ^  V 


where 


A  =  (  Un+l  )/5  L  n  P  (2n+i  /2)  }  ] 


-.2n*l)/4  *  2r.»i  ;  6  *  VQ]’ 


q  =  -I  /8  -  :  1  -'64  *  1  /?',  ^ 


B  =■  (q-n^/^p-ql  A 


C  ■••  (p-n)/(p-q)  A 


Conbi r.at ion3  cf  n  and  P  that  satisfy  the  relation 
n  »  -1/8  ^  (i/6i»  ♦  l/(2P))^''^ 

:a'ose  the  solution  for  A  in  Equation  3-‘i9  to  become  indefinite.  For  t 
special  oases,  particular  solutions  to  Equation  must  be  calculated, 

example,  with  n  =  1,5  and  P  «  ^/2^,  the  solution  is 

V  -  -(M7/100)  [X^'^^ln  X  -  X^''^/5] 

Next,  the  more  general  case  without  the  small  angle  assumption  is  cor.sid 
General  case:  wave  angle  not  necessarily  small.  On  gently  si  op  in 
Peaches,  the  wave  angle  at  breaking  is  usually  small  due  to  wave  refract 
cut  this  is  not  always  the  case.  Liu  and  Dalrymple  (1973)  and  Kraus  and 
Sasaxi  ( 1979a}  show  that  the  breaking  wave  angle  has  a  significant  effec 
tne  magnitude  of  the  longshore  current  and  the  shape  of  the  current 
iist'-icution.  The  method  of  Kraus  and  Sasaxi  is  followed  to  include  tne 
effect  cf  wave  angles  on  the  longshore  current  distriouti on. 
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3y  Snell’3  law  and  Equation  3-16 


*.  V  V V. 

s:. 


3in0  -  C/C^3in8^ 
b  b 


1  /2 

=  (h/h.  )  3in6^ 

b  b 


Using  a  trigonometric  identity,  the  cosine  of  the  wave  angle  may  be  written  as 

co3e  =  :i  -  (b/n,,)  sin^e,.)’'^  :-5s 

b  b 


i*  'i>  "j 


Equations  3-55  and  3-06  express  the  sine  and  cosine  of  the  local  wave  angle 
terms  of  water  depth  and  constants. 

Using  Equation  3"56,  the  local  wave  stress  in  the  surf  zone,  Equation 
3-20,  becomes 


•/*/*■  V/ 


Ty  >  tanB*l/8  p  g^''^Y^/(h^)  sine./C 


b  b 


S/3h  (i-(f^/)T^)3in^e^)'  ''2s 


*  ■-  *»  •Ji 

•  V  V  ■■ 
V  S  •• 


'aking  the  derivative  and  simplifying.  Equation  3”57  becomes 


'  =  p  g  Y“  tan8*/i6  sine 

D  D 


-.n* '  '  (1 


ih/h^Jsin  0^) 


-  (h/h^)sin^9j^/(l-(h/h^)3in^9^)  ' 


The  expression  for  the  local  wave  stress  outside  the  surf  zone  remains  tr.e 
same.  The  lateral  mixing  stress  was  unaffected  by  the  small  angle  assampt: 
therefore  Equation  3~25  is  still  valid.  The  bottom  friction  stress  :Equat; 
3-314)  becomes 


irp  V 


%  g’^^h^'^"’'^^V(h^)^""'^i*(h/h.  }  sin^e,  ) 


g^ H/h^'"^  (’♦ih/h^)  sin^e^) 


X  V  X. 


X  >  X, 


/•/■/’iT] 

*.  *•  ••  ^ 
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The  stress  balance  inside  the  breaKer  line  from  Equations  3-5' 
and  3-34  is 


:>  » 


-p  g  Y  tans  /16  3in0^  h 

b 


(2n-l/2)/^^  ^ (2n-3/2) 
"b 


:  (  4n*l  )  ; :  -  'h/h,  )  sin^e, 

0  D 


',/h,  ;  3;n^9^/;'-( h/h,  ;  sin^S^)'"'^:  = 
0  b  D  b 


1/2  ( n- 1 t  r  1  /7 1 

or/2  Y  g  '''*^13'  ’d/dx  'x  h“  “'9v/3x) 


’  /  o  (  n- 1  /2 )  n- 1  )  0 

'nag  -  V  Y  ■'^V(h^)"^  ^l*X3in-9j 


Equation  3"60  is  nondimensionalized  and  simplified,  resulting  in 


P  3/3X(X^'^*^^^^3V/3X)  -  V  X^""’ ^^\l*X3in‘0,)  =» 

D 

-:4n-’;/5  X'^‘'’'’'^^:(l-Xsin^9,,}’'^^  -  X/(Un*l)  sin^d^/ ( 1 -Xsir.^e^  } ' ' 

b  b  b  ■ 


'axing  the  derivative  in  Equation  3~6l  yields 


P  (n*3/2)  X  V  ♦  P  X^V  ’  -  V  (1*X3in'9^',  =  '3-52 

-:4n*i}/5  X‘^[(1-X3ln^9^)’ -  X/(4n+l)  3in^9^/(  l -Xsin'^g,, )  ' ''^: 

0  b  b 

The  quantity  in  square  brackets  in  Equation  3-52  is  approximated  by  a  :i'',or.ia'. 
expansion  truncated  to  second  order 


[1  -  ( 4n+3)X ( 2(  4n*l  ) )  X  sin'^S.- 

b 


(  4n*5  )/(8(  4n+1  )  )  X^  3in"*9,] 
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Krau3  and  Sasaki  (1979a,  I979b)  obtained  a  solution  to  Equation  3-62  in  tne 
form  of  an  infinite  series  by  retaining  all  orders  of  the  binomial  expansion 
of  the  breaking  wave  angle.  Truncation  past  second  order  is  here  considered 
to  be  sufficiently  accurate  and  allows  a  more  convenient  solution  for 
engineering  application.  Equation  3-62  then  becomes 


P  (n*3/2)  X  V  *  P  X^V"  -  V  (1  *  X  sin^9,) 

b 


,%  .'iVv 


-b,  X  »  b^  X-  *  b^  X- 


b,  «  (iln-l  )/5 


b,  >  ('*n*3)  3in^9^/(i0) 
2  b 


r.jr^TT; 


b,  -  («n*5)  sin  0^/(‘JO) 
3  t> 


s  ,“».>■ 


Equation  3~b4  is  a  second-order  nonhomogeneous  differential  equation.  The 
solution  to  Equation  3“6^  for  the  region  shoreward  of  the  breaker  line  is 
approximated  by  a  power  series  truncated  to  second  order 


•  •  V.*  ■ 

:-r\.  * 

.■■rN*.  •> 


A  X"  >  A  x"*’  ♦  A  X^*^ 
o  1  2 


3  X^  ♦  3  X^*’  *  3.  x'^*“ 

o  1  2 


where 


Aq  -  [(4n-1  )/5]/[1-(n',(2n*0.5)P] 


A^  •  [-( Un*3)3in^9^/1 0  -  A^sin^  9  ] /1 1  - ( n* 1 i { 2n*l . 5 ) ? ! 


A2  -  [-(‘In>5)3in^9^/‘i0  -  A,3in^  9  ] /[  1  - (  n*2)  (  2n-2 . 5 i  P ; 


3^  -  (Bg3in^9j^)/:(p*1)(p*n*l  .5)P-1  : 


'Iv'^  v'..-  •i'liv'lvvxv’ju 


--**  v*  *■ 


•.v.i 
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■pe  stress  balance  seaward  of  tne  PreaKer  line  ^Epjaticn  3--iM  with  tr.e  wave 


reignt  approximated  by  Green's  Law  Equation  3'^-'  nondimensi  oral  form  ts 


a 


i  G  -X  X"'"  3V  3x:  -  x''"^-;  -  -  X  sm^e  ' 


ilo. Gating  tr.e  oe'-ivative  in  Equation  3-06  and  rearranging  tr.e  te.^-ms  yields 


a  x“  v’  -  5/^  a  X  V’  -  Vi’  ♦  X  sin'e^;,  -  o 


Equation  3-62  is  a  homogeneous  second-order  differential  equation.  The 
solution  of  Equation  3-6'’  is  approximated  by  a  power  series  truncated  to 


second  order 


x*’  ♦  x'’"’*  x*^ 


q  =  -(1/8)  -  (l/eii^l/Q) 


C,  =■  3in^9^/[(q-1 )  ;q-3X‘*)Q-'' 


C2  “  C.|  3ln“9(j/[(q-2)  (q-7/4)a-'’ ] 


Expressions  for  the  coefficients  3  from  Equation  3‘oo  am 


■or.  rq^at; 


:-o3  are  obtained  by  equating  the  current  and  the  derivative  of  the  ourr-^ 


.nside  and  outside  the  surf  zone  at  the  breaKer  line.  The  general 


ihe  longshore  current  distribution  truncated  to  second  orde'-  in  tne 


ireauirg  wave  angle  is 
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wnere  is  the  wind  stress  in  the  y-direction,  tr.e  drag  coefficient 

Barrett  197"'),  w  is  the  wind  speed,  and  $  is  the  incident  wind  directicn 
Wilson  I960).  Birkemeier  and  Dalrymple  present  a  nearshore 

circulation  model  that  includes  the  effect  of  wind  stress.  The  analytical 
solution  of  the  stress  balance  with  the  addition  of  the  wind  3t"es3  as  giver, 
in  Equation  3~70  with  lateral  mixing  neglected  is 

V  =  sin’e^)  l"'Un*-i  '/5  X^\'  )  ^ -Xsin'e  )  - 

J  D 

Xsin^a.  /  ,  1 -Xsin'e,  '  '  “  ) 

o  6  -  ' 

?  '  n- ; ' 2 ) 

*  T  C..W  3initi/(v  Y  .gh)^  c.)  1/X' 

J  0  D  I 

Tne  solution  becomes  indefinite  near  tne  shoreline.  This  problem  could  oe 
overcome  by  representing  the  fluid  flow  and  wind  stress  in  the  swash  zone  ne'e 
accurately.  Such  a  task  is  beyond  the  scope  of  this  report. 

Discussion  of  Results 

The  main  points  discussed  in  this  section  are:  a'  the  effect  of  the 
«ave  neight  decay  power  law  on  the  longshore  current  distribution,  b'  a 
ccmpariscn  of  the  longshore  current  model  with  data,  and  o'  the  mode. 

1 .-itations. 

The  longshore  current  .model  does  not  reduce  exactly  to  tne  model  of 
Long-et'Kiggins  ) 1970a,  1970b)  for  n  =  i  and  small  incident  wave  angle 
because  the  form  of  the  lateral  viscosity  coefficient  fellows  Madsen  et  al. 

't^S)  instead  of  Longuet-Hlggins,  This  is  not  a  f u.nda.mental  difference,  an: 
.t  will  not  be  considered  in  the  discussion. 

The  effect  of  the  exponent,  n  ,  on  the  longshore  current  profile  is 
shewn  in  Figure  3-2.  Increasing  n-vaiues  steepen  the  current  p''ofile  and 
inc''ea3e  the  distance  from  the  shoreline  to  the  maxim-un  velocity.  T.nerefore, 
I':-'  mild  teach  slopes  and  large  values  of  X  the  n-value  will  be  la-ge,  tne 
:  .r’"ent  distribution  will  be  more  peaked,  and  the  l:cati:n  of  tne  maxi-..un 
:  .r'"ent  -ill  be  closer  to  the  breaker  line.  Typical  values  of  n  "ange  f'':- 
t;  2.1.  Figures  3-3  and  3~-‘  show  tne  effect  of  varying  tne  ?-val^e.  i 
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small  value  of  P  indicates  the  bottom  friction  stress  dominates  tr.e  late”: 
mixing  stress,  and  a  large  value  of  =  indicates  the  lateral  mixing  stress 
dominates  the  bottom  friction  stress.  Larger  values  of  P  flatten  and 
broaden  the  current  profile.  The  value  of  n  =  l  in  Figure  3-3  corresponds  •; 
t^.e  Longuet-Higgins  model.  For  r.  =  1  and  P  =  0,  the  current  profile  is 
t”! angular  (note  the  discontinuity  at  the  breaker  line  as  explained  earlier 
f  1”  tne  no-*ateral  mixing  case,.  In  Figure  for  n  =  '.5  and  =  =  I,  tr.e 

d”df.l:»  .3  concave  upward  with  tne  same  discontinuity  at  the  bresKer  line. 


.g..”e3 


■3,  and  3"^  are  for  an  approximately  zero  incident  wave  angle, 


.-.ply.ng  tne  higher  order  wave  angle  effect  was  omitted. 

Tne  effects  associated  with  increased  wave  angle  are  shown  in  Figure  3 
Tne  nondimensional  current  decreases  with  increasing  breaking  wave  angle.  Ti 
valie  of  V  in  Figure  3'5  for  an  incident  breaking  wave  angle  of  30  degrees 
n  =  ’.0  and  ?  -  0.05)  is  30  percent  lower  than  for  an  incident  breaking  wav' 
a'gle  of  0  degrees.  Also,  the  location  of  the  maxim-urn  current  is  closer  to 
tne  shoreline  witn  increasing  breaking  wave  angle. 

Tne  longshore  current  model  developed  herein  was  compared  to  laborator 
!sta  r”om  yizugucni  et  al.  (*9"'3)'  and  to  the  model  of  Longuet-Higgins.  A 


the  laboratory  data  is  found  in  Kraus  and  Sasaxi 


As  n;t' 


:y  •(''=..5  and  Sasaki,  the  position  of  the  maximum  current  velocity  va.-ies 
: :  ns  -  der  abl  y ,  and  is  therefore  a  good  parameter  for  correlating  tne  ;ur'’ent 
model  prediction  with  the  laboratory  observations.  The  position  of  tne 
-a:<imu.m  velocity  is  used  to  determine  P  ,  Cr  ,  and  T  ,  given  '•  ,  end  r, 
Tne  met.nod  employed  by  Kraus  and  Sasaki  to  esti.T.ate  the  para.meter5  f.-or.  tne 
cate  IS  used  with  the  additional  step  to  determining  n  from  ''  and  ta-; 

Figures  3-6,  3-7,  3-8,  and  3-9  illustrate  the  fit  of  tne  iongs-.c.-e 
Current  model  and  the  Longuet-Higgins  model  to  the  Minz'-gucni  et  ai.  cat  a. 
Tne  current  velocity  is  normalized  by  the  maximum  velocity.  Inside  f-.e  s.-: 
zone,  the  model  fits  the  data  well.  Near  the  shoreline,  tne  effect  of  t-e 
cower  law  decay  can  be  seen  in  the  slightly  concave  upward  snape  of  one 
:”ofile.  "Taple  3-1  gives  the  values  of  and  0^  oalouiated  oy  fito.-.g  O'- 
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and 


values 


the  model  are  s.igntlv 


gne 


n  •  1. 00 
P  =  0.05 


(no  0b  correction) 
0b-  10.0* 

0b-  20.0* 
^  0b-  30.0* 
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'ig-jre  3-5.  Longshore  current  distribution  shoving  the 
iependence  on  the  breaking  wave  angle  =  j.Z-  'no  h: 

:rder  angle  effects),  10. OO,  20.0'^,  30.3^T  n  =  I.IC,; 

=  =  0.05) 
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Fig'j.'-e  3-6.  Comparison  of  the  longshore  current  distribution  from  Case 
experimental  results  of  Mizuguchi  et  al.  (1978),  the  present  model,  and 
the  model  of  Longuet-Higgi.ns  (1970b) 
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Figure  3-'’.  Comparison  of  the  longshore  current  distribution  from.  Case 
experiment.3l  results  of  Mizuguchi  et  al.  (1973),  the  present  model,  ano 
tne  model  of  Longuet-H iggi ns  ■1970b) 
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Figure  3“8.  Comparison  of  tne  longshore  current  distribution  from  Case  3 
experimental  results  of  Mizuguchi  et  al,  ;i978),  the  present  model,  and 
tne  .model  of  Longuet-Higgins  i’970b) 
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'igure  3"9.  Comparison  of  the  longshore  current  distributitn  from  Case 
experimental  results  of  Mizuguchi  et  al.  (i978),  the  present  mcdel,  a-.d 
ine  model  of  Lor.guet-H iggins  (1970&) 
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outside  the  surf  zone  slightly  better  than  the  Longuet-Kiggins  .-nodel,  but  the 
departure  from  the  data  is  still  large.  The  reasons  for  tnis  difference  are 
not  known,  but  may  be  due,  in  part,  to  the  accuracy  of  measuring  the  breaker 
positioh,  the  accuracy  of  measuring  low  current  velocities,  and  the  effect  of 
circulation  in  the  enclosed  wave  basin. 

The  data  set  of  Mizuguchi  et  al.  does  not  rigorously  test  tr.e  current 
hcdel.  The  data  were  collected  on  a  slope  of  V'O.iJ,  so  the  expected  values 
;f  r.  are  close  to  l.O  as  shown  in  Chapter  II.  An  n- value  of  ^  .  0  ''educes  tne 
bro.<en  wave  height  to  a  linear  function  of  the  water  depth,  and  the  current 
".cdel  '-educes  to  a  truncated  version  of  the  model  of  Kraus  and  Sasawi,  or  to 


Table  3-I 

Mizuguchi  et  al.  Longshore  Current  Data 


Case  1 


(cm/s 


idicates  the  vaiue  for  the  Longuet -Higgins  model 


cf  Longuet-Higgins  if  the  incident  wave  angle  is  small.  Beacres 
■d  States  consisting  of  C.2-mm  sand  typically  have  slopes  in  tne  r 


of  1/140  to  1/70.  The  n-valuea  for  these  milder  slopes  would  be  greater  than 
1.0,  and  the  effect  of  the  power  law  wave  height  decay  on  the  current  profile 
would  be  more  pronounced. 

Table  3“''  illustrates  a  more  subtle  point.  Earlier  work,  corresponding 
to  a  value  of  n  »  1 ,  may  require  somewhat  different  values  of  Cr  and  P  to 
fit  the  data.  Again,  the  data  of  y.izugucni  et  al.  does  not  test  this  point 
rigorously  because  the  values  of  n  are  close  to  unity. 

The  application  of  the  longshore  current  model  is  limited  net  only  by 

the  assumptions  listed  in  Table  i-2,  but  also  by  the  truncation  of  the  power 

series  solution  to  second  order.  The  effect  of  the  truncation  increases  as 

tne  values  of  9^  ,  ?  ,  and  n  increase.  The  effect  of  varying  these 

parameters  over  typical  ranges  ;s  examined.  The  value  of  6  is  limited  to 

less  than  approximately  30  degrees  because  of  the  truncation  of  the  bottom 

friction  stress.  Typical  values  of  P  range  from  0.01  to  C.IO,  and  typical 

values  of  n  range  from  1.0  to  2.0.  For  a  value  of  ?  equal  to  0.5  and 

5^  •  30'^,  the  maxim'ura  difference  between  the  infinite  power  series  solution 

and  the  truncated  solution  is  only  6  percent  and  at  6  -  20  the  difference 

D 

reduces  to  less  than  '  percent.  Figure  3~5  shows  the  truncated  solution  for  n 
=  ’.0  and  P  »  0.05,  and  Figure  3"10  shows  the  infinite  series  solution.  For  a 


?-value  of  O.’O,  the  difference  between  the  infinite  series  and  truncated 
series  solutions  is  11  percent  for  9^  •  30'^  and  ^  percent  for  0,^  -  20"  . 
an  r.-value  of  2.0  (P  *  0.05)  the  difference  between  the  infinite  series  and 
tr’ur.cated  series  solutions  are  37  percent  for  9^  *  30"  and  u  percent 
for  9^  -  20^  .  Figure  3~'I1  shows  the  truncated  series  solution  and  -igure 
"2  snows  the  infinite  series  solution  for  n  •  2.0  and  ?  •  0.25.  In  summary 


estimates  the  infinite  series  well  for  relatively  small  values  of 
ci.t  caution  should  be  used  applying  the  model  for  large  values  of 
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CHAPTER  IV;  WAVE  SETUP 


Wave  setup  and  setdown  are  the  change  in  the  ^ean  water  level  due  t« 
excess  momentum  in  the  x-direction.  In  the  surf  zcne  there  is  normally  a 
setup  of  the  water  level,  whereas  seaward  of  the  breaker  line  there  is  a 
setdown.  The  wave  setup  in  tne  longshore  current  model  is  approximated  c 
altering  the  beach  slope  as  given  by  Equation  This  chapter  describ 

the  derivation  of  the  wave  setup  from  the  equation  of  motion  in  the  x- 
direction  (Equation  3-IJ)  based  on  the  power  law  description  of  the  broken 
height.  Although  this  form  of  the  wave  setup  is  not  included  in  the  long 
current  mode*,  it  is  an  application  of  the  power  law  wave  heignt  decay. 

The  equation  of  motion  in  the  x-direction  becomes 

C  =  og(h»n)c(n/3x*3S  i3x 

with  the  only  x-directed  stress  being  the  principle  component  of  the  '•adi 
stress,  .  The  quantity  n  is  the  time-mean  water  surface  elevation 

to  wave- induced  momentum.  The  mean  flux  of  momentum  across  a  plane  x  » 
constant  is 

S  -  3/2  E 

XX 

Shallow  water  ( Longuet-H iggins  and  Stewart  ’964).  Substituting  in  tne 
energy  density  given  by  Equation  3-''0,  Equation  4-2  becomes 

S  -  3/'!6  p  g  H'" 

XX 

and  tne  momentum  balance  (Equation  )  expands  to 

j  -  p  g  (h  ♦  ni  3n/3x  ♦  3/3x  ^2/^6  p  g  H" 


Shoreward  of  the  breaker  lire,  the  power  law  wave  height  decay  i: 
applied  to  describe  the  broken  wave  height,  and  Equation  a-d  becomea 


p  g  (h  ♦  n)  3n/3x  ♦ 


3 '  lo  P  g  3/3x  (Y^;h^  ';h  ♦  'i;/;h^  ♦ 


iliulatlhg  the  derivative  in  Equation  y-5  and  simplifying,  yields 


1  =  ,  h  ♦  ''1  3  n  /  3  X 


0  -  f  2n-2 '  -  >'  2n-  i  ) 

B-'io  "'"''^h,  ♦  n.  )  '2n)  'h  ♦  n;  i3n/3x  +  3h/3x) 

b  b 


Eearr  ar.gi  r.g  the  terms  gives 


3n/3x  »  /1 1 


♦  K/(h  ♦  n) ^ :  3h/3x 


-  '  P  p  -  ✓ 

K  ■  3 '  h  *  r  ) ' 

b  w 


'I'.g  for  r  sy  treating  h  ♦  n)  as  a  single  variable  ar 


n,  -  -(h  ♦  n)  ■  'K(2n-'  : 


-'.ere  C  is  a  constant  of  integration. 

Seaward  of  the  breaker  line  the  energy  flux  Is  constant  ar 
.3  given  by 


r  »  -H  '  0  r. ' 


•  V  •  V  V  V  >V  V  V  V.V .  V.W.  V.  V.  iT,  rf-  -  -  J-  J-  -  -f  • 


/  vi  -'.vy  .-I’.* ' 


In  shallow  water  (Longuet-Higgins  and  Stewart  196U).  This  is  referred  ■:  o  as 
"setdown"  because  it  is  a  depression  of  the  mean  water  surface.  Equating  tne 
solutions  for  n  seaward  and  shoreward  of  the  breaker  line  at  the  breaker 


line  and  noting  ^  yields  the  solution  for  the  integration  constant 

C  >  (Y“h^/'!6}r'-3n>^/3)  »  ^n  ♦  1  ]/ {2n  -  i  ) 

D 

The  solution  of  Equation  u-3  tecones 
-  ^  2  n  -  ■'  ) 

r  -  "  n}'  /C.K(2n*  ^)) 

'Y  n.^/'t':.'-3n^‘^/3;  *  ^n  ♦  'y[2n  -  M 

The  solution  for  n  must  oe  found  iteratively  because  Equation  is 

implicit.  For  the  special  case  of  n  *  1  (linear  wave  height  decay,,  n  can 
te  expressed  explicitly 

r  -  (3y'^/.:3  ♦  3Y^)):-n  *  hj^/i6(-Y^  ♦  UO/3}'/  --•2 

F:n  this  special  case  the  setup  is  a  linear  function  cf  the  water  cectn. 
values  of  h  greater  than  one,  the  profile  of  the  wave  setup  is  concave 
townward . 

"igures  U-’  and  4-2  illustrate  the  effect  of  the  exponent  f-'or.  tne  to-e'' 
ia-  -ave  neight  decay,  h  ,  on  the  wave  setup  profile.  Figure  -■'.z-a  tne 

p-ofi.e  cf  the  wave  setup  calculated  from  Equation  ( r.  =  an: 

Equation  4-12  (n  -  1.00)  for  the  small-scale  experimental  nun  of  Stive 
■335'.  The  setup  measured  by  Stive  is  also  plotted.  Figure  --2  snows 
o-ofile  of  the  wave  setup  calculated  from  Equation  '  n  -  '.oS  ant 
Equation  a-i2  (n  -  i.OC)  for  the  large-scale  experimental  r.n  of  Stive. 

Again,  the  setup  measured  by  Stive  is  also  plotted.  The  n-values  ^seq 
Equation  4-12  were  calculated  with  Equation  2-3  from  the  teach  slope  ' 
ant  the  measured  ^-values. 


The  setup  profiles  calculated  usir.g  both  Equation  4-n  and  EquatiDn 
overestimate  the  wave  setup,  but  the  calculated  setup  based  on  the  power  law 
wave  height  decay  represents  the  data  better  than  the  calculated  setup  based 
cn  linear  wave  height  decay.  The  difference  between  the  two  calculated 
profiles  is  greatest  at  the  point  of  maximum  setup,  which  is  the  critical 
point  in  most  engineering  studies. 

In  a’antmary,  the  wave  setup  is  calculated  based  on  the  power  law 
expression  of  the  broxe.n  wave  height  developed  in  Chapter  II.  For  the  liritec 
amount  of  setup  data  exam.ir.ed,  the  calculated  setup  based  on  the  power  law 
-ave  neignt  decay  describes  the  trend  of  the  measurements  better  than  t^e 
oaloulated  setup  cased  on  linear  wave  height  decay.  5oth  express. o.ns 
: ve’' es ti mate  the  data. 
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CHAPTER  V:  CONCLUSIONS 


The  p'jrpose  of  this  report  was  to  develop  an  analytical  model  of  the 
longshore  current  based  on  a  power  law  expression  for  the  broken  wave  heigr.t 
in  the  surf  zone.  The  model  was  intended  to  be  an  improvement  over  present 
models  based  on  linear  wave  height  decay.  For  use  as  an  engineering  tool,  the 
model  was  to  be  as  general  as  possible,  including  the  effect  of  wave  setup, 
finite  wave  angle,  and  lateral  mixing. 

An  empirical  power  law  expression  for  the  broken  wave  height  was 
developed  based  on  seven  independent  data  sets  consisting  of  135  experimental 
runs.  The  exponent  of  the  power  law  expression  is  a  function  of  the  beach 
slope  and  the  ratio  of  wave  height  to  water  depth  at  wave  breaking.  From  the 
data,  typical  values  of  the  exponent  range  from  1.0  to  2.0.  High  values  of 

the  exponent  correspond  to  mild  beach  slopes,  small  ratios  of  the  wave  height 

to  water  depth  at  wave  breaking,  and  concave  upward  wave  height  profiles.  For 
an  exponent  equal  to  1.0,  the  broken  wave  height  reduces  to  a  linear  function 
of  the  water  depth.  In  previous  longshore  current  models,  a  linear  wave 
height  decay  was  assumed  for  all  beach  slopes  and  breaker  height  to  breaker 
depth  ratios.  The  power  law  decay  is  shown  to  represent  the  wave  height  decay 
profiles  significantly  better  than  linear  decay.  The  power  law  decay 
expression  also  compared  favorably  to  the  more  complex  decay  model  of  Tally  et 
al.  '1985a,  1985ti),  To  use  the  power  law  expression  in  a  predictive  mode,  the 
ratio  of  wave  height  to  water  depth  at  wave  breaking  must  be  estimated.  This 

'’atio  is  best  estimated  by  the  expression  of  Singamsetti  and  Wind  i'98C)  or 

the  expression  of  Sunamura  (I98l), 

The  longshore  current  model  is  based  on  the  momentum  balance  in  the 
longshore  direction.  Many  simplifying  assumptions  are  made  in  the  model  in 
order  to  provide  a  solution  in  a  form  for  practical  use.  The  driving  force  of 
the  longshore  current  is  the  local  wave  stress  which  is  calculated  using  the 
::ncept  of  radiation  stress.  The  lateral  mixing  stress,  caused  cy  horizontal 
•  .''0  jlent  eddies,  redistributes  momentum.  Flow  of  water  along  the  bo' tom  is 
-  . ]‘ed  0/  the  bottom  friction  stress.  The  derivation  of  the  longshore 
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current  follows  the  method  used  by  Longuet-Higgins  (1970a,  ’970b).  The  effect 
of  incident  wave  angles  is  included  in  the  form  presented  by  Kraus  and  Sasaki 
v’9'’9a,  1979b).  Wave  setup  is  accounted  for  by  altering  the  beach  slope.  The 
longshore  current  is  expressed  as  a  power  series  in  the  wave  angle  truncated 
tc  second  order.  Wind  stress  is  not  included  in  the  general  solution, 
although  some  examination  was  made  of  its  effect. 

The  *ongshore  current  model  was  compared  to  laboratory  data  from 
'^izuguchi  et  al .  '19'^3(  and  the  model  of  Longuet-Higgins  (1970b).  The 
longshore  current  model  represents  the  data  well,  although  it  appears  to 
underestimate  the  current  speed  seaward  of  the  breaker  line.  There  is  some 
doubt  about  the  validity  of  the  data,  however,  for  the  seaward  region.  The 
longshore  current  model  follows  the  trends  of  the  data  slightly  better  than 
tne  Longuet-Higgins  model,  but  the  data  set  is  not  a  rigorous  test  of  the 
model.  The  experiment  was  performed  on  a  steep  beach  slope,  so  the  expected 
exponent  in  the  power  law  decay  expression  is  close  to  1.0,  reducing  the  wave 
height  decay  to  approximately  a  linear  function  of  water  depth.  The  range  of 
incident  wave  angles  for  which  the  model  can  be  applied  is  limited  by  the 
truncation  of  the  power  series  solution,  but  covers  the  useful  range  of 
realistic  breaking  wave  angles. 

The  mean  wave  setup  and  setdown  are  derived  from  the  momentum  balance  in 
the  shore  normal  direction  based  on  the  power  law  wave  height  decay.  The 
profile  of  the  wave  setup  is  concave  downward  for  an  exponent  in  the  power  law 
wave  height  decay  greater  than  unity,  whereas  the  setup  calculated  from  linear 
wave  height  decay  is  linear.  The  estimated  setup  based  on  the  power  law  wave 
height  decay  represents  the  setup  data  collected  by  Stive  (1985)  better  than 
the  estimated  setup  based  on  linear  wave  height  decay.  Both  calculated 
estimates  of  setup  overestimate  the  measurements. 

The  understanding  of  wave  height  decay  and  longshore  currents  gained 
from  this  investigation  suggests  areas  of  future  study:  a)  collection  of 
additional  longshore  current  data  to  verify  the  longshore  current  model  over’ 
the  range  of  typical  beach  slopes,  wave  height  to  water  depth  ratios  at  wave 
breaking,  and  incident  wave  angles;  b)  collection  of  additional  longshore 
our^’ent  data  to  quantify  the  friction  coefficient  and  the  eddy  viscosity,  so 
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the  model  could  be  better  applied  in  a  predictive  mode:  c)  application  of  t 
longshore  current  model  to  predict  the  distribution  of  sediment  transport 
across  the  surf  zone;  d)  application  of  the  power  law  wave  height  decay 
directly  to  other  wave  energy  problems  in  the  surf  zone  (e.g.,  sediment 
transport  and  wave  setup);  and  e)  extension  to  random  wave  breaking. 


REFERENCES 


3asco,  D.  R.  1982.  "Surf  Zone  Currents,"  Miscellaneous  Report  No.  82-7  (I), 
C.  S.  Army,  Corps  of  Engineers  Coastal  Engineering  Research  Center,  Fort 
Belvoir,  VA. 

sasco,  D.  R.,  and  Coleman,  R.  A.  1982,  "Surf  Zone  Currents,"  Miscellaneous 
Report  No.  32-7  (II),  U.  S.  Army,  Corps  of  Engineers  Coastal  Engineering 
Research  Center,  Fort  Belvoir,  VA. 

Rattjes,  J.  A.  1975.  "Surf  Similarity,"  Proceedings  of  the  ^^th  Coastal 
Engineering  Conference,  American  Society  of  Civil  Engineers, 
pp  •466  -‘J30. 

Battjes,  J.  A.,  and  Janssen,  J.  P.  F.  M,  1979.  "Energy  Loss  and  Setup  due  to 
Breaking  of  Random  Waves,"  Proceedings  of  the  ^6th  Coastal  Engineering 
Conference ,  American  Society  of  Civil  Engineers,  pp  569-587. 

Birkemeier,  W.  A.,  and  Dalrymple,  R.  A.  1975.  "Nearshore  Water  Circulations 
Induced  by  Wind  and  Waves,"  Symposium  on  Modeling  Techniques,  American  Society 
of  Civil  Engineers,  pp  1062-1081. 

Bowen,  A.  J.  1969.  "The  Generation  of  Longshore  Currents  on  a  Plane  Beach," 
Journal  of  Marine  Research,  Vol  27,  No.  1,  pp  206-215. 

Bowen,  A.  J,,  Inman,  D.  L. ,  and  Simmons,  V.  P.  1968.  "Wave  'Set-Down'  and 
Setup,"  Journal  of  Geophysical  Research,  Vol  73,  No.  8,  pp  2569-2577. 


Collins,  J.  I.,  and  Wier,  W.  1969.  "Probabilities  of  Wave  Characteristics  in 
the  Surf  Zone,"  Tetra  Tech  Report  #TC-1U9. 

Dally,  W.  R.,  Dean,  R.  G.,  and  Dalrymple,  R.  A.  1985a.  "A  Model  for  Breaker 
Decay  on  Beaches,"  Proceedings  of  the  19th  Coastal  Engineering  Conference, 
American  Society  of  Civil  Engineers,  pp  82-98. 

_ .  1985b.  "Wave  Height  Variation  Across  Beaches  of  Arbitrary 

Profile,"  Journal  of  Geophysical  Research,  Vol  90,  No.  C6,  pp  Il9i7-li92h. 

Dean,  R.  G.  1977.  "Equilibrium  Beach  Profiles:  U.S.  Atlantic  and  Gulf 
Coasts,"  Ocean  Engineering  Report  12,  Department  of  Civil  Engineering, 
University  of  Delaware,  Newark,  DE. 

Dean,  R.  G.,  and  Dalrymple,  R.  A.  198U.  Water  Wave  Mechanics  for  Engineers 
and  Scientist,  1st  ed.,  Prentice-Hall,  Inc.,  Englewood  Cliffs,  NJ. 

DivoKy,  D.,  Le  Mehaute,  B.,  and  Lin,  A.  1970.  "Breaking  Waves  on  Gentle 
Slopes,"  Journal  of  Geophysical  Research,  Vol  75,  No.  9,  pp  1681-1692. 


N  ^  -  J 


»■  a 


.  V  V- 


>  >r 

N.  N.  u" 


%r  v.“, “-w  •«  *.  ^  % 


Galvin,  C.  J.  1967.  "Longshore  Current  Velocity:  A  Review  of  Tr.eory  ar.t 
Gata,"  Review  of  Geophysics,  Vol  5,  No.  3,  pp  297-3CA. 

_ _ .  1969.  "Breaker  Travel  and  Choice  of  Design  Wave  Height," 

Journal  of  the  Waterways  and  Harbors  Division,  American  Society  of  Civ:l 
Engineers,  Vol  95,  No.  'WWZ,  pp  175-200. 

Garrett,  J.  R.  1977.  "Review  of  Drag  Coefficients  over  Oceans  and 
Continents,"  Monthly  Weather  Review,  Vol  105,  pp  915-929. 

Goda,  Y.  1970.  "A  Synthesis  of  Breaker  Indices,"  Transactions  of  the  Japan 
Society  of  Civil  Engineers,  Vol  2,  Part  2,  pp  227-230. 

Harris,  T.  E.  W.  1969.  "Nearshore  Circulations;  Field  Observiti ons  and 
Experimental  Investigations  of  an  Underlying  Cause  in  Wave  Tanks,"  Symposium 
on  Coastal  Engineering,  Stellenbosch,  South  Africa. 

Horikawa,  K.,  and  Kuo,  C.  1967.  "A  Study  on  Wave  Transformation  Inside  the 
Surf  Zone,"  Proceedings  of  the  Uth  Coastal  Engineer ing  Conference,  American 


Society  of  Civil  Engineers,  pp  217-233. 

Hwang,  L.,  and  Divoky,  0.  1971.  "Breaking  Wave  Setup  and  Decay  on  Gentle 

Slope,"  Proceedings  of  the  i 2th  Coastal  Engineering  Conference,  American 


Society  of  Civil  Engineers,  pp  377-389. 

James,  I.  D.  197^.  "A  Non-linear  Theory  of  Longshore  Currents,"  Estuarine 
and  Coastal  Marine  Science,  Vol  2,  No.  3.  PP  235-2U9. 

Jonsson,  I.  G.  1967.  "Wave  Boundary  Layers  and  Friction  Factors," 

Proceedings  of  the  iQth  Coastal  Engineering  Conference,  American  Society  of 
Civil  Engineers,  pp  127-198. 

Jonsson,  I.  G.,  Skovgaard,  0.,  and  Jacobsen,  T.  S.  1975.  "Computation  of 
Longshore  Currents,"  Proceedings  of  the  I9th  Coastal  Engineering  Conference, 
American  Society  of  Civil  Engineers,  pp  699-719. 

Keely,  J.  R.,  and  Bowen,  A.  J.  1977.  "Longshore  Variation  in  Longshore 
Currents,"  Canadian  Journal  of  Earth  Science,  Vol  19,  No.  8,  pp  l897-’905. 

Kraus,  N.  C.,  and  Sasaki,  T.  0.  1979a.  "Influence  of  Wave  Angle  and  Lateral 

Mixing  on  the  Longshore  Current,"  Marine  Science  Communications,  Vol  5,  No.  2, 

pp  91-126. 

_ .  1979b.  "Effects  of  Wave  Angle  and  Lateral  Mixing  on  the 

Longshore  Current,"  Coastal  Engineering  in  Japan,  Japan  Society  of  Civil 
Engineers,  Vol  22,  pp  59-79, 

Kuo,  C.  T.  1965.  "A  Study  on  Wave  Transformation  After  Breaking," 
unpublished  M.  S.  Thesis,  Department  of  Civil  Engineering,  University  of 
Tokyo,  Tokyo,  Japan.  Un  Japanese) 


Le  Mehaute,  B.  1963.  "On  the  Nonsaturated  Breaker  Theory  and  the  Wave  Run 
Up,"  Proceedings  of  the  8th  Coastal  Engineering  Conference,  American  Society 
of  Civil  Engineers,  pp  77-92. 

uiu,  P.  L.-F.,  and  Dalrymple,  R.  1978.  "Bottom  Frictional  Stresses  and 
Longshore  Currents  Due  to  Waves  with  Large  Angles  of  Incidence,"  Jo-urnal  of 
'Marine  Resea'"ch,  Vol  36,  No.  2,  pp  357-375. 

Longuet-Higgins,  M.  S.  l97Ca.  "Longshore  Currents  Generated  Py  OPliquely 
Incident  Sea  Waves,  i,"  Journal  of  Geophysical  Researcn,  Vol  75,  No.  33,  op 

_ .  ''9’^Ob.  "Longshore  Currents  Generated  by  Obliquely  Incident  Sea 

Waves,  2,"  Journal  of  Geophysical  Research,  Vol  75,  No.  33,  PP  6"90-p8C1 . 

Longuet-Higgins,  M.  S. ,  and  Stewart,  R.  W.  1962.  "Radiation  Stress  and  Mass 
T-’ansport  in  Gravity  Waves,  with  Application  to  'Surf  Beats',"  Journal  of 
Fluid  Mechanics,  Vol  13,  No.  U,  pp  1481-504. 

_ .  1963.  "A  Note  on  Wave  Setup,"  Journal  of  Marine  Research,  Vol 

:• .  No.  ’ ,  pp  14-10. 

_ .  '964.  "Radiation  Stress  in  Water  Waves;  A  Physical  Discussion 

w.th  Applications,"  Deep  Sea  Research,  Vol  ll,  No,  4,  pp  529-562. 

'^adsen,  G.  S.,  Cstendorf,  D,  W.,  and  Reyman,  A.  S.  1978.  "A  Longshore 
Durrent  Model,"  Symposium  on  Technical,  Ehvironmental ,  Socioeconomic,  and 
Regulatory  Aspect  of  Coastal  Zone  Management  Coastal  Cone  ’"3,  American 
Scciety  of  Civil  Engineers,  pp  2332-2341  . 

Maruyama,  K.,  Sakakiyama,  T,,  Kajima,  R.,  Saito,  S. ,  and  Shimizu,  T.  '983. 
"Experimental  Study  on  Wave  Height  and  Water  Particle  Velocity  Near  the  Surf 
Cone  Using  a  Large  Wave  Flume,"  Civil  Engineering  Laooratory  Report  No. 

The  Central  Research  Institute  of  Electric  Power  Industry,  Chiba, 
Japan.  'in  Japanese) 

McCcwan,  J,  1891,  "On  the  Solitary  Wave,"  Philosophical  Magazine,  Journal 
Science,  London,  Edinburgh,  Dublin,  Vol  32,  No.  5. 

McDougal,  W.  G.,  and  Hudspeth,  R.  T.  1986.  "Influence  of  Lateral  Mixing  on 
Longshore  Currents,"  Ocean  Engineering,  Vol  13,  No.  5,  pp  419-433. 


Meadows,  G.  A.  1977.  "Time  Dependent  Fluctuations  in  Longshore  Currents," 
^;"oceeding3  of  the  15th  Coastal  Engineering  Conference,  American  Society  of 
Civil  Engineers,  pp  660-680. 

Melville,  W.  K.  1977.  "Wind  Stress  and  Roughness  Length  Over  Breaking 
Waves,"  Journal  of  Physical  Oceanography,  Vol  7,  pp  702-710. 


A' 


■'-ew 

_  • 

k.  j-  . 


k.  is’-v 


V!iller,  I.,  and  Freund,  J.  E.  1977.  Probability  and  Statistics  for 
Engineers ,  3<3  ed.,  Prentice-Hall,  Inc.,  Englewood  Cliffs,  N  J . 

Mizuguchi,  M.,  Oshima,  Y.,  and  Horikawa,  K.  1978,  "Longshore  Currents," 
P-oceedings  of  the  25th  Japanese  Conference  on  Coastal  Engineering,  Japan 
Society  of  Civil  Engineers.  (in  Japanese) 

Miz'iguchi,  M.  I98l.  "An  Heuristic  Model  of  Wave  Height  Distribution  in  Surf 
Cone,"  Proceedings  of  the  1'^th  Coastal  Engineering  Conference,  American 
Society  of  Civil  Engineers,  pp  278-289. 

Nakamura,  M.,  Shirashi,  H.,  and  Sasaki,  Y.  1967.  "Wave  Decaying  Due  to 
Breaking,"  Proceedings  of  the  lOtn  Coastal  Engineering  Conference,  American 
Society  of  Civil  Engineers,  pp  25^-253. 

Saeki,  H.,  and  Sasaki,  M.  1973.  "A  Study  of  the  Deformation  of  Waves  After 
Breaking  (1),"  Proceedings  of  the  20th  Japanese  Conference  on  Coastal 
Engineering,  Japan  Society  of  Civil  Engineers,  pp  559-564.  (in  Japanese) 

Sasaki,  M.,  and  Saeki,  H.  1974.  "A  Study  of  the  Deformation  of  Waves  After 
Breaking  (2),"  Proceedings  of  the  21st  Japanese  Conference  on  Coastal 
Engineering,  Japan  Society  of  Civil  Engineers,  pp  39-44.  (in  Japanese) 

Sawaragi,  T.,  and  Iwata,  K.  1975.  "On  Wave  Deformation  After  Breaking," 
P-cceedings  of  the  I4th  Coastal  Engineering  Conference,  American  Society  of 
Civil  Engineering,  pp  481-499. 

Singamsetti,  S.  R.,  and  Wind,  H.  G.  1980.  "Characteristics  of  Breaking  and 
Shoaling  Periodic  Waves  Normally  Incident  on  to  Plane  Beaches  of  Constant 
Slope,"  Report  M1371,  Delft  Hydraulics  Laboratory. 

Stive,  M.  J.  F.  1985.  "A  Scale  Comparison  of  Wave  Breaking  on  a  Beach," 
Coastal  Engineering,  Vol  9,  pp  151-158. 

Street,  R.  L.,  and  Camfield,  F.  E.  1967.  "Observation  and  Experiments  on 
Solitary  Wave  Deformation,"  Proceedings  of  the  10th  Coastal  Engineering 
Conference,  American  Society  of  Civil  Engineers,  pp  234-293. 


yy.'-'j.'* 
.  • 


v'.V- 

s^:v‘ 


0  •  . 


Suhayda,  J.,  and  Pettigrew,  N.  R.  1977,  "Observations  of  Wave  Height  and 
Wave  Celerity  In  the  Surf  Zone,"  Journal  of  Geophysical  Research,  Vol  82,  No. 
9,  pp  1  419-1  424. 

Sunamura,  T.  1981.  "A  Laboratory  Study  of  Offshore  Transport  of  Sediment  and 
a  Model  for  Eroding  Beaches,"  Proceedings  of  the  I7th  Coastal  Engineering 
Conference,  American  Society  of  Civil  Engineers,  pp  1051-1070. 

Svendsen,  I.  A.  1984.  "Wave  Heights  and  Setup  in  a  Surf  Zone,"  Coastal 
Engineer-. ng,  Vol  8,  pp  303“329. 


V, 


-f.  "-.d 
'.V,  V. 


36 


.-.w'.-'ii* 

'•■nW 


,  ■  V-  .■ 


.hornton,  E.  B.  1971.  "Variation  of  Longshore  Current  Across  tne  Surf  Cone,' 
Proceedings  of  the  I2th  Coastal  Engineering  Conference,  ATierioan  Society  of 
Civil  Engineers,  pp  291-308. 

Thornton,  E.  B.  ,  and  Guza,  R.  T.  1983.  "Transf  ormation  of  Wave  Heigr.t 
Cistribution,"  Journal  of  Geophysical  Research,  Vol  38,  pp  5925-5938. 

Van  Dorn,  W.  G.  1977.  "Setup  and  Run-up  in  Shoaling  Breakers,"  Proceedings 
of  the  'Sth  Coastal  Engineering  Conference,  American  Society  of  Civil 
Engineers,  pp  733-751. 


Weggel,  J.  R.  1972.  "Maximum  3reat<er  Height,"  Journal  of  tne  Waterways, 
Harbors,  and  Coastal  Engineering  Division,  American  Society  cf  Civil 
Engineers,  Vol  93,  No.  W'W9,  pp  529-5^8. 

Wilson,  3.  W.  i960.  "Note  on  Surface  'Wind  Stress  over  Water  at  Low  and  High 
Wind  Speeds,"  Jo'urnal  of  Geophysical  Research,  Vol  65,  No.  10, 

PP  337"-3332. 


^ 

I  M 


APPENDIX  A:  WAVE  HEIGHT  DECAY  DATA 


S  =  slope 

T  =  wave  period  (s) 

HO  =  deepwater  wave  height 
DB  =•  water  depth  at  wave  breaking 
H  =  wave  height 
D  =  water  depth 

:ii)  =  heights  and  depths  in  entimeters 
T.)  =  heights  and  depths  in  meters 
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.71 

5.37 

.65 

5.42 

.59 

4.97 

.53 

4.52 

.47 

4.07 

.41 

3.61 

.35 

1 /30.0 

T  = 

2.0 

H 

D/DB 

15.50 

1  .00 

15.00 

.95 

1  4.06 

.90 

13.60 

.85 

11.25 

.80 

9.3"^ 

.75 

7.50 

.70 

■'.03 

.65 

6.66 

.60 

6.10 

.55 

5.16 

.50 

4.69 

.45 

4.22 

.40 
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3.28 
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1 /65.0 

T  - 

2.0 

H 

D/DB 

24.00 
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.91 

13.10 

.89 

’4.70 

.84 

1  3.50 

.73 

3.71 

.73 

■’.2’ 

.67 

I'. 21 

.  6l 

0.82 

.56 

6.31 
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.41; 

15.22  DB  >  21.3  (cm) 


16.33  DB  =  25.0  (cm) 


19.26  DB  -  27. 'I  (cm) 
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5.90  DB  =  10.9  (cm) 
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3.91 
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3.91 
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.65 
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.77 

6.88 
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H 
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1 .00 

’1.26 

.93 

9.40 

.87 

^.31 

.65 

6.44 
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7.34 
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.81 
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.75 

6.41 

.69 

7.49 

.63 

5.66 

.56 

5.37 

.49 

4.02 

.41 

3.82 

.35 

3.66 

.27 

3.19 
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i  a 
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H 

D/DB 

6 . 36 

1  .00 
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4.61 

.73 

4.32 

.67 
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.38 

2.  31 
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1  .  4 

H 

D/DB 
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1.00 
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.37 

4.75 

.73 

3.83 

.56 
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.42 

2.37 

.27 
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.13 

■'20.0 
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H 

D/DB 

-  _  -5 

1.00 

6.22 

.71 

3.93 

.46 

2.31 

.21 

/20.0 

T  - 

1  .4 

H 

D/DB 

9.42 

1 .00 

9.00 

.96 

6.93 

.83 

c  _  O'' 

.67 

^.22 
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2.69 
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1/20.0  T 
H 

10.26 

9.47 

6.97 

5.68 

4,41 

4.51 

3.56 

2.95 

’ /20.0  T 
H 

n  .  1  ij 
9.06 
5.38 
5.37 
3.99 
3.49 
3.25 
2.87 

1/20.0  T 
H 

11.63 
9.55 
3.29 
5.54 
-*  .  43 


'/2C.0  T 
H 


1.4  HO 
D/  DB 
1.00 
.73 
.69 
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.44 
.36 
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.16 

1  .4  HO 
D/DB 
1.00 
.31 
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T 
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1 1  .98 

.80 

8,03 
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.27 
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H 
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5.00 
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2.80 
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T  >  2.2  HO  -  9.46 

DB  -  10.5 

( cm) 

H 

D/DB 

12.34 

1.00 

11.39 

.88 

9.00 

.76 

6.32 

.64 

5.93 

.50 

4 .  a6 

.38 

3.71 

.26 

2.69 

.1  4 

1  '20.0 

T  -  2.2  HO  -  9.83 

DB  .  12.3 

( cm) 

H 

D/DB 

'  ’  .  32 

1.00 

’  0.’6 

.86 

1  0.62 

.76 

7.99 

.65 

6.22 

.54 

5. ’3 

.43 

4.52 

.33 

3.82 

.22 

2.73 

.12 

o 

o 

f\J 

T  -  2.2  HO  -  10.12 

DB  •  ’2.9 

■  cm ; 

H 

D/DB 

’  ’  .  60 

’  .00 

•  O.-O 

.88 

3  •  0  0 

.69 

6.^2 

.59 

6 .  26 

.-9 
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5.4 

4.2- 

3.2 


39 

29 


1/20.0 

T  - 

2.2 

H 

D/DB 

11.69 

.98 

1  2.81 

1  .00 

11.39 

.91 

1  ’  .06 

.31 

10.  30 

.74 

9.71 

.72 

".30 

.65 

6.54 

.56 

5.33 

.46 

5,07 

.37 

1^.50 

.28 

3.33 

.19 

1  /20.0 

T  » 

2.2 

H 

D/DB 

1  2.  31 

1 .00 

1  1  .25 

.91 

9.29 

.81 

8.27 

.72 

6.72 

.62 

5.73 

.47 

5.34 

.42 

5.53 

.40 

'■i .  34 

.30 

3.56 

.21 

2.59 

.13 

'  '20.0 

2.2 

H 

D/DB 

'  2. 

1  .00 

■  ;.o5 

.92 

’2.20 

.34 

9.98 

.74 

9.29 

.66 

7.37 

.57 

5.9-" 

.49 

5.77 

.41 

5.98 

.3^ 

^^.95 

.33 

^.52 

.26 

3 , 4  *5 

.’5 

2.4'^ 

.03 

10.57  DB  =  13.5  (cm) 


1.34  DB  «  13.3  (cm) 


12.14  DB 


1/20.0  T  - 

2.3 

H 

D/DB 

15.58 

1 .00 

14.94 

.92 

1  4.79 

.85 

11.00 

.74 

10.85 

.67 

8.19 

.59 

".''0 

.52 

6.-2 

.44 

6.  ’  3 

.36 

5.27 

.29 

“^.52 

.24 

3.-2 

.  1  7 

2.81 

.09 

’/50.0  T  > 

1.3 

H 

D/DB 

'  0.60 

1.00 

’  0.34 

.93 

’0.23 

.96 

10. ’2 

.94 

9.65 

.92 

9.01 

.39 

3.35 

.37 

3.37 

.35 

3.06 

.83 

.30 

7,26 

.30 

■:.68 

i 

.76 

5.33 

."5 

- .  25 

.74 

5.5’ 

5 .  ’  4 

.67 

- .  82 

.64 

4 .  ’  3 

.63 

-.45 

.62 

4.13 

.59 

4.1  3 

.58 

4.35 

.57 

4.03 

.55 

4.03 

.53 

3.98 

.52 

- .  92 

.50 

4.03 

.  49 

3 . 60 

.43 

• .  5  ’ 

.46 

3.32 

.  -4 

3.55 

.  ..  1 

13.65  DB  -  16.5  (cm) 


=  1 0.  31  DB  =■  16.4  (cm) 
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2.47 
2.47 
^.55 
6.50 
1  0.53 
1  0.27 

1  3.00 


/25.O 

H 

.  30 
.’9 
.  94 
.94 
.77 
.59 


4.3  HO 
D/DB 
.20 
.34 

.50 
.  6u 
.-^3 

.85 

1  .00 

’  .6  HO 
0/DB 
.07 
.15 
.20 
.25 


=  1 3.0c  DB  =  13.3 


=  ’5. 40  DB 


m 

m 

“w  % 


.s 

*•  ,*‘.N 

v.-’.v 


2.4  HO 
D/DB 
.10 
.19 
.32 
.43 
.51 
.59 
.64 


1 4,40  DB  =  16.9  (cm) 


ro  Lf>  LT-  O 


s  >  1/25.0  T  -  3,4  HO 

H  D/DB 

2.48  .15 

2.48  .28 

4.37  .40 

5.43  .51 

5.90  .79 

9.09  .91 

•'.30  1.00 

S  =  •  25.0  T  =  4.3  HO 

H  C/DB 

3.2'  .14 

3.43  .23 

''.:5  .40 

- .  76  .50 

5.59  .65 

5.04  .-3 

9. ■"5  .91 

' ’ . 90  1.00 

2  =  '  '  ’ 2 . 0  T  =  ’ . 6  HO 

H  D/  DB 

2.03  .10 

6.24  .27 

5.77  .35 

6.40  .51 

11,08  .70 

'4.32  .84 

’5.60  1.00 

3  =  ’/1 2.0  T  =  2.4  HO 

H  D/DB 

5.46  .17 

5.21  .46 

7.24  .60 

12. •’’0  .88 

12.70  1.00 

S  =  1/12.0  T  "  3.4  HO 

H  D/DB 

.70  .20 

.77  .37 

.^7  .53 

.95  .69 

.47  ,88 

'4.30  1.00 


1  1 .80  DB  =  1  1  .1  (o.n) 


1  1 . 90  DB  =  1  1  . 1  ( cm ) 


15.60  DB  =  18.3  (cm) 


12,70  DB  =  10.8  (cm) 


1  4.80  DB  =  9.3  (cm) 
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APPENDIX  B:  LONGSHORE  CURRENT  COMPUTER  PROGRAM  AND  SAMPLE  RUN 

LONGSHORE  CURRENT  BASED  ON  POWER  LAW  WAVE  HEIGHT  DECAY 

This  program  calculates  the  longshore  current  based  on  an  empirical  power 
expression  for  the  wave  height  decay  in  the  surf  zone, 

H  =  r  *hb»(h/hb)»»n 

where : 


H 

i  3 

the 

h 

is 

the 

hb 

is 

the 

Hb 

is 

the 

p 

is 

the 

n 

is 

the 

The 


exponent  (typical  range  ’.0  to  2.0) 

The  exponent  n  is  a  function  of  the  beach  slope  and  the  breaker  index 
exponent  may  be  input  directly  or  calculated  in  the  program  from  the  beach 
slope  and  breaker  index.  Other  inputs  include  the  parameter  P  expressing 
the  relative  importance  of  lateral  mixing  and  bottom  friction  (typical  range 
1.01  to  0.10)  and  the  wave  angle  at  breaKing  (typical  range  0.0  to  30.0 
degrees ) . 


Do  you  want  the  program  to  calculate  the  power  law  exponent?  (Y  or  N' 

y 

Input  the  beach  slope,  breaker  index  (e.g.  0.02,0.3) 

1.02, 3. 78 

Input  the  parameter  P 
3 . 05 


(nput  the  breaking  wave  angle  (degrees) 


=  1 .7U  p  =  0.05  Breaking  wave  angle  =  10.0 
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N  N 


ongshore  current  Distribution 

V  is  the  dimensionless  longshore  ourrent  speed 
X  13  the  dimensionless  distance  from  tr.e  mean  sho'-ei 
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Analytical  Model  of  the  Longshore  Current  Based  on  < 
Power  Law  Wave  Height  Decay  Including  the  Effect  of  ‘ 
Large  Incident  Wave  Angles 

JANE  SMITH  29  OCT  86 

lefir.  itions  of  variables  and  arrays 

N  -  power  law  exponent,  input  or  calculated  internally 
from  the  beach  slope  and  breaker  index 
?  -  parameter  expressing  the  relative  importance  of 
lateral  mixing  and  bottom  f-'iction  (0.0  1  to  O.'C) 
THETA  -  breaking  wave  angle  in  degrees 
M  -  order  of  the  solution  {2'. 

THETAR  -  breaking  wave  angle  in  radians 
SINTB2  -  sine  of  the  breaKing  angle  squared 
From  Equation  3"^9: 


p  D  - 

p 

,  GO  -  q  .  A(  1  ) 

-  Ao ,  A( 2 '  -  A 

’  ,  a: 

3.1  -  A  2 

BETA 

V  1  . 

)  -  Bo,  BET  A (2) 

-  31 ,  beta: 3) 

-  B2 

SUMA 

SNMA 

A  (  • 

I)  -  Co,  DELTA (2)  -  Cl,  DELTA 
SA,  SUMS  -  SB,  SUMC  -  SC 

SA  '  ,  SNM3  -  SB  '  ,  SNMC  -  SC 

(3)  - 

C2 

3(1) 

- 

bl  ,  B(2)  -  b2. 

c :  1  )  -  d  ,  C  ( 

2)  - 

c2 

X ; :  1 

■ 

no  n  d i me  hs i ona 1 
x  /  X  b 

distance  from 

the 

shoreline 

V : :  ) 

nondimensional 
V  /  vo 

longshore  cur 

rent 

speed  , 
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DIMENSION 
DIMENSION 
DATA  I  N  0  /  ' 
REAL  N 
M»2 

V  M  -  ^  4  1 

TYPE  60 
F0RMAT( //  , 

*  1 X  ,  ' HEIGHT 
TYPE  61 
FORMAT ( / 1 X 
based ' , 

*  '  on  an 
»  ,  ’  height 
•//’x.'wher 

*  '  9  X  ,  '  h  b  i 
* ’ X , ' he i ght 

*  '  9  X , ' n  is 
TYPE  62 


A(25),3ETA(25),DELTA(25), 3(25  1,0(25) 
X( 1 00) , V; 100) 

N ' / , I  YES/ ' Y ' / 


1  X  LONGSHORE  CURRENT  BASED  ON  POWER  LAW 
DECAY ' ) 

,  'This  program  calculates  the  longshore  : 


empirical  power  law  expression  for  the 
decay  in  the  surf  z one ,  '  ,  /  /  5 x  ,  ’ H  «  gamm 
e:  H  is  the  wave  height',/9x,'h  is  t 

s  the  water  depth  at  break! ng',/9x,'Hb 
at  breaKing’,/9x, 'gamma  is  the  breaKer 
the  exponent  (typical  range  i.O  to  2.0 


he  wave’ 
amma*h.o* 


e  water  dec' 
is  t  n  e  w  a  V  e 
index  H  C  n : 


^  O  •  T- 1 


V,  V, 


s"  s'^V  V 
j-  V 

’KK-^ 

•  M  *  »  • 


* 


FORMATS'  The  exponent  n  i 3  a  function  of  the  beach  slope  a 
•/'  the  breaker  index.  The  exponent  may  be  input  directly  cr 
*/'  calculated  in  the  program  from  the  beach  slope  and  brea^e 
*/'  index.  Other  inputs  include  the  parameter  P  expressing 
♦/'  the  relative  importance  of  lateral  mixing  and  bottom', 

*/'  friction  (typical  range  O.Oi  to  O.iO)  and  the  wave  angle' 
»/'  at  breaking  (typical  range  0.0  to  30.0  degrees).') 

TYPE  63 

F0RMAT(/'/,lX,'Do  you  want  the  program  to  calculate  the  power 
•■x.'law  exponent?  (Y  or  N)') 

READ  6y.:AN3 
F  IR^AT { A ’  ) 

:f(ians.eq.:ng)  g:  to  to 

TYPE  66 

FI P. MAT(1X, 'Input  the  beach  slope,  breaker  index  'e.g.  0.02,1 

READ  •.SDCPE.BINDEX 

N=0.:32-0.0096/SLOPE*0.657*3INDEX*O.OP3*BINDEX/SLOPE 
GO  TO 
TYPE  72 

FORMAT (/' 1 X  ,' I  nput  power  law  exponent,  n') 

READ  »  ,  N 
TYPE  76 

FORMAT (/ 1 X  ,' I  nput  the  parameter  P') 

READ  *  ,  P 
TYPE  77 

FCRMAT(/1X, 'Input  the  breaking  wave  angle  'degrees)') 

READ  •, THETA 
TYPE  78,N,P,THETA 

FORMAT)//, IX,  'n  -  ',FU.2,2X,'P  =  ',Fa,2,2X,'3reaKing  wave', 
*^x, 'angle  •  '.FR.I,///) 

F  A  :  N  •  1  .0 

fa:n2*i .0 

THETAR>3.16159*THETA/130.0 

3INT32»(SIN(THETAR))»»2 

^°=(-(2.0*N*1  .0)/4.0)*3QRT(((2.0*N*1  .C)/4.0)*»2*1  .0/P) 

;Q-(-1  ./8.)-SQRT(1  .0/64. 0*1  .0/?) 

A(1}.(4.0»N  +  l.0)/5.0/(l.0-P»N*(2.0*N*0.5'') 

BETA (  1  )  -  1  .  0 
DELTA (  1  )  -  1  .  0 
SUMA- A ( 1 ) 

3UMB-BETA ( 1  ) 

SUMC-DELTA( 1 ) 

SNMA-N* A( 1 ) 

SNMB-PP»BETA( 1 ) 

3NMC-QQ*DELTA ( 1 ) 


00  10  I  -  1  ,  M 
I  I  -  I  ♦  1 
Z  -FLOAT (  I  ) 

fa:n-facn*z 

F  A  C  N  2  -  :  2  .  1  ♦  Z  -  .  1)  *  F  A  C  N  2 


"L- '^-  •-  ».- 


IF(FACN2.LT.l  .0)FACN2-1  .0 

A(II)-(-(2.0*Z+4.0*N*1)*FACN2)/:5.0«?.0»»Z*FACN) 

A(II)-A{II)*SINT92*»2-SINT32*A(II-1) 

A(i:)-A(II)/(1.0-P*(N-Z)*(2.0»N*2*0.5)) 

BETA(II)-SINTB2»BETA(i:-n 

beta(i:)-beta:ii)/(?»(p?  +  z)*(p?*-z-n>o.5;-'!.;) 


:*0.25)-1  .01 


OELTA ' II ) »SINT32»DEZTAv  :i-1  ) 

D  E  LT  A  (  I  I  )  -  D  E  LT  A  ;  1 1  )  /  ;  ?  •  f  -  Z  1  »  '  QQ  - 

SUyA-S’JV!A*A(  II  ) 

SUyB  =  S'JM3*BETA  :  II  ) 

S'jyc  =  SL:yc*oELTA(  ii  i 
SNyA  =  SNyA  +  (  N-^Z  )  *  A  (  I  I  ) 
3Ny3=SNyS*:PP+Z) »3ETA(II) 
SNyC=SNMC*(QQ-Z)*DELTA(II) 

CONTINUE 


BM),(SNyA*3UMC-SUMA»SNMC)/(3UMB*SNyC-SNy3*SUMC) 
c;  1  )  =  ';3NMA»3UM3-SUMA*SNMB)/(SUMB*3NyC-SNyB»SUyC) 

00  20  1=2, MM 
3 ; I }  =  3 ! 1 ) *BETA :  I  ) 

C  (  I  )  =C  (  1  )  *DELTA ( I  ) 

CONTINUE 

OC  1 00  J= 1  ,  50 

X( j; =FL0AT( J-1  )/50. 

XU*50)=X(J)+1  .0 
V ;  J :  =  0 . 0 
:  j  -  5  0 )  =  0 . 0 
00  101  K=1 ,MM 

:  =  FL0AT ( K ) -  1  . 0 

V:J)=V(J)*A(K)*X(J)**(N-Z)*3(K)*X(J)*»(PP*Z; 

V(J*50)-V(J*50)+C(K)»(X(J+50))*»(QQ-Z) 

CONTINUE 
CONTINUE 
TYPE  43 

FORMAT ( 1 X Longshore  Current  D i s t r i & u t  i  o n '  1 
TYPE  114 

F0RMAT(5X,'V  Is  the  dimensionless  longshore  current  speet' 
*/5x,'X  Is  the  dimensionless  distance  from  the  mean  shore  1: 
DO  45  1-1,10 
K  1  =  (  I  -  1  )  *  1  0  ♦  1 
K2-K 1 -9 

TY°E  50,(V(K),K-K1,K2) 

TYPE  51  ,  (X(K)  ,K-K1  ,K2) 

TYPE  52 
CONTINUE 

FORMAT (  1 X  V  -  '  ,  ’ 0F7  .  3  ) 

F  0  R  y  A  T  :  i  X  ,  '  X  =  1  0  F  7  .  3  ) 
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APPENDIX  D:  NOTATION 

aA'-< 

• 

\ 

\  <3y> 

V  ^ 

lateral  mixing  stress 

average  bottom  friction  stress 

.•  if,  W,i 

1 

V 

1 

wave  celerity 

• 

! 

\ 

drag  coefficient 

h  -'’ 

f^'iction  coefficient 

:  '? 

group  wave  celerity 

s‘  •. 

vy  :■'! 

1 

•  S  V 

1 

wave  energy  density  per  unit  surface  area 

• 

1  f  ^ 

Doriolis  parameter 

'  "  X 

energy  flux  in  the  onshore-direction  per  unit  distance 

1 

parallel  to  shore 

^  3 

gravitational  acceleration 

• 

'  H 

wave  height 

•fc 

water  depth 

V^v  f 

•“’•  K 

wave  neight  at  wave  breaking 

'  ‘v!’ 

• 

;  -b 

water  depth  at  wave  breaking 

U 

"0 

deepwater  wave  height 

;  < 

decay  coefficient  in  Dally  et  al.  (1985a,  1985b)  wave  neight 
decay  model 

i 

n 

-■  'o 

/v'v'.- 

* 

deepwater  wavelength 

beach  slope 

• 

V  N  \ 

I  r. 

exponent  in  power  law  wave  height  decay  expression 

'  ? 

parameter  expressing  the  relative  importance  of 

k 

lateral  mixing  and  bottom  friction  (including  the 

•v 

» 

effect  of  wave  setup)  in  the  longshore  current  model 

• 

1 

r  - 

parameter  expressing  the  relative  importance  of 

w*. 

.V  .S 

lateral  mixing  and  bottom  friction  (excluding  the 

n\v',n 

effect  of  wave  setup)  in  the  longshore  current  model 

w-y-v**/! 

'  XX 

lilCuti  k  lU  V^Ui  ^  1 1 U  dlU  dVJf  d  ^Xdi*^ 

principle  component  of  radiation  stress 

k 

• 

1 

: 

mean  flux  of  y-momentum  parallel  to  the  shore  across  a 

-  - . 

plane  x  -  constant,  component  of  radiation  stress 

1/_  ✓ 

w’ 

> 

wave  period 

^ 

time 

• 

1 

.'-'.'V*.'-  ."5 

•»-0 

•  >’>.■ 

•■•Vv'-.-l 

129 

“\.V 

.  P. 

. 

onahore  current  component 
■wave  orbital  velocity 

nondimensional  longshore  current  speed,  v/v^ 
longshore  current  component 

maximum  longshore  current  speed  for  special  case  of 
n  =  I  and  lateral  mixing  stress  neglected 

dimensionless  onshore  coordinate, 

onshore  coordinate 

location  of  'wave  breaking 

alongshore  coordinate 

wind  speed 

angle  of  bottom  'with  the  horizontal 
constant  in  the  expression  for  e, 

U 

ratio  of  wave  height  to  water  depth  at  wave  breaking 
energy  dissipation  rate 
lateral  viscosity  coefficient 
wave  setup 

wave  setup  at  wave  breaking 

incident  wave  angle 

incident  wave  angle  at  wave  breaking 

the  constant  -rr 

density  of  water 

incident  wind  angle 

local  wave  stress 

wind  stress 


